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1.0 INTRODUCTIONS 


1.1 CTS THERMAL ANOMALIES 

The Communication Technology Satellite (CTS) was launched into an 
equatorial geosynchrouous orbit in January, 1976 and is stationed at 116° 
longitude. After completing nearly four years of operation, the CTS mis- 
sion has recently been terminated. 

The major payload on the CTS is the Transmitter Experiment Package 
(TEP) which consists of a high power travelling wave tube (TWT), a power 
processing system (PPS) and a variable conductance heat pipe system (VCHPS) 
shown schematically in Figure 1-1. Heat produced in the tube collector is 
directly radiated to space, whereas power dissipated in the PPS is first 
conducted to the TEP saddle and the spacecraft south panel from where it 
is radiatively rejected. 

The temperature control of the TWT is provided primarily by a variable 
conductance heat pipe (VCHP)/Radiator system which transports and rejects 
most of the heat dissipated in the tube body. Mounting details of the tube 
body to heat pipe evaporator saddle and the spacecraft south panel are 
sketched in Figure 1-2. As shown there are several interfaces through which 
the heat must flow to reach the heat pipes. Normally, heat reaching the 
aluminum baseplate is mostly absorbed by the heat pipes and is transported 
axially to the VCHP radiator. A fraction of it, however, is conducted to 
the spacecraft south panel where it is radiated to space. 

The VCHP system, shown in Figure 1-3, consists of three stainless 
steel /methanol dual artery heat pipes whose variable conductance is achieved 
with the help of external cold wicked reservoirs loaded with a 10 percent 
helium/90 percent. nitrogen gas mixture. Figure 1-3 shows the six positions 
in the VCHP/radiator system at which temperatures are measured in flight. 

The Transmitter Experiment Package is instrumented with four additional 
temperature sensors, two located on the Multistage Depressed Collector (MDC) 
and two sensors located on and near the tube body. 

The CTS is the first spacecraft to rely on heat pipes to control the 
thermal performance of a major onboard system. The Transmitter Experiment 
Package has performed satisfactorily except on four occasions in 1977, 

March 16, March 23, April 11, and September 10 when measured tube body 
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temperatures displayed sudden, rapid increases not normal for the TEP 
operating conditions and inconsistent with the design of the thermal con- 
trol systems. The anomalous behavior of the TEP in these four occasions 
are referred to as the "CTS Thermal Anomalies” of days 75, 82, 101, and 
253 respectively. 

During all the anomalies, the continuous increase of the tube body 
temperatures was reversed by reducing the output power without damage to 
TEP. Operation of the TWT at the reduced power resulted subsequently in 
the recovery of the thermal control system with no evidence of degradation 
or changes in its performance. 

After the fourth anomaly (day 253), the TEP was successfully operated 
at somewhat reduced power (below the open artery capacity of the heat pipes), 
and no further anomalies were detected. 

1.2 REVIEW OF PREVIOUS WORK 

Following the first occurrence (day 75), the CTS anomalies were the 
subject of investigations both at NASA Lewis Research Center (LeRC) and at 
TRW in an effort to understand the cause of the thermal failures, focussing 
particular attention on the VCHPS since it is in the thermal path from the 
tube body to space. The work performed at LeRC and preliminary studies 
conducted at TRW are documented in References 1 and 2, respectively. 

LeRC's investigation included the use of flight-type TEP components 
in ambient and vacuum ground tests, analytical studies, flight data review, 
and CTS on-orbit tests. The objectives were to determine the most probable 
cause of the anomalies and to identify procedures or TEP operating modes 
which would preclude damage to the TWT in the event of recurring anomalies 
or of continuous degraded capacity of the thermal control system. 

From the review of flight data for both normal and anomalous days the 
following observations about the anomalies were made: 

1) All occurred after relatively long periods at constant, 
although differing, TWT RF power levels or heat rejection 
rates from the tube baseplate. 

2) All the anomalies took place near the Spring and Fall equinox 
periods during which sun angle with respect to the VCHPS 
radiator is relatively small and during the orbit quadrant 
where the CTS is almost directly between Earth and Sun in 
which case the VCHPS surface becomes increasingly shadowed by 
the spacecraft itself. 
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3) All the anomalies, except on day 253, were characterized by being 
preceded by temperatures measured at the extremity of heat pipe 
number 1 (HPl with sensor HPT5. See Figure 1-3) being at or 
below -98 C, the methanol freezing point. 

4) All the thermal anomalies were preceded or accompanied by a 
sudden abnormal but small increase in the difference in 
measured temperatures between the adiabatic sections of heat 
pipe number 3 and 1, aTj_ 3 . This increase in aTi _3 was 
followed by a second increase, and at the outset of unstable 
tube body temperature rise, aT 3 _i observed to drop to a 
lower level. 

Tests were performed with a TWT almost identical to the one outboard 
CTS. The results indicated that the anomalies could not have been caused 
by thermal interface failure inasmuch that: a) it is improbable that a 
failed interface could recover its integrity completely after an anomaly, 
and b) it was demonstrated that an interface failure results in tube body 
temperature rates higher than those observed during the anomalies. 

Vacuum tests that most closely approximated anomalous temperature 
profiles were those in which the tube base plate cooling was abruptly 
reduced from that required for thermal equilibrium at normal operating 
temperatures. These test results suggested that the reduction in heat 
transfer capability of the variable conductance heat pipe system is, to 
a high decree of certainty, the cause of the thermal anomalies. 

Room-ambient tests at LeRC with a flight- type VCHPS indicated that 
aT 3 _i increases are caused by depriming of heat pipe 1 (HPl) or heat pipe 
(HPl) and heat pipe 2 (HP2). The temperature difference between the adia- 
batic sections of heat pipe 3 (HP3) and HPl was found to depend on the 
priming states of the heat pipe arteries. The measured temperature dif- 
ferences aT 3 _i were: 

a) aT 2_2 = I.IC, all heat pipes are primed 

b) ATgj = 5. 1C, HPl is deprimed 

c) ATg_^ = 7.8C, HPl and HP2 deprimed 

d) aT 2_2 = 2.3C, all heat pipes deprimed 
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These values for aT 3 _i are consistent with the temperature levels observed 
during anomalous days and tend to confirm the postulate that the anomalies 
are initiated by depriming of HPl, followed by dspHming of HP2, and the 
onset of tube body unstable temperature rise caiised finally by depriming 
of HP3. 

On day 253, while LeRC personnel were performing experiments on-board 
CTS, the fourth anomaly occurred. This event gave the unique opportunity 
to establish the stable VCHPS heat rejection capacity beyond which the 
tube body temperature becomes unstable. It was found the open artery 
capacity of the VCHPS is approximately 106 watts. 

The comprehensive investigative work performed at LeRC made signifi- 
cant contributions to the understanding of the CTS thermal anomalies. The 
most significant contribution is the fact that convincing evidence was 
made available to establish that the thermal anomalies were caused by 
unexpected reductions in the VCHPS heat transfer capability resulting from 
depriming of the heat pipe arteries. 

The first preliminary investigation of the CTS thermal anomalies con- 
ducted at TRW is documented in Reference 2. This investigation, performed 
in direct support of LeRC efforts, focussed particular attention at the 
identification and analysis of potential artery depriming mechanisms. 

Selection of deprimiJig mechanisms was based on the premise that they 
must be consistent with three key observations: 

1) All the anomalies occurred during the eclipse season. 

2) These anomalies are sporadic. Similar conditions on successive 
days can yield an anomaly one day but not the next. 

3) The anomalies are triggered by depriming of arteries in HPl. 

In the early study, among a number of postulated mechanisms, only 
three were identified which could be associated with eclipse seasons: 

I) Condenser freezing 

II) Marangoni flow 

III) Gas evolution in the arteries due to rapid chilling 
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It was argued that under certain conditions freezing of CIS arteries could 
lead to depriming. Arterial failure could result from the transfer of mass 
from the evaporator to the frozen condenser where it would freeze and not 
be available for circ;ilation. This transfer of mass could be effected by 
vapor diffusion, by liquid pumping in fillets due to reduction of specific 
volume as it cools and/or by Marangoni flow. The fact that the CTS-type 
heat pipes will not fail due to condenser freezing when operated under no 
load or at high heat loads, led to the argument that the condenser freez- 
ing mechanism will work if the pipe is operated at the right load, not too 
high, not too low. However, several factors implied that, although con- 
denser freezing can lead to artery depriming, it could not be the primary 
cause of the anomalies. 

These factors are the fact that on day 253 an anomaly occurred before 
the radiator portion of HPl approached freezing conditions, and the fact 
that with all pipes primed HP3 should freeze and deprime first since it 
carries least load and its radiator portion radiates from both sides, how- 
ever, HPl failed first. An additional factor is the fact that after the 
anomalies the VCHPS recovers before the radiator begins to warm up. It 
was thus argued that if sufficient excess fluid was frozen to cause artery 
depriming there would be none available to rewet the evaporator. There- 
fore, condenser freezing was ruled out. 

Marangoni Flow (surface tension flow due to temperature gradients 
along the a VCHP gas-blocked region) could induce a flow along fillets and 
excess fluid reservoirs toward the condenser end along the vapor-liquid 
interfaces. No flow, however could be induced in the wick and arteries 
which have structure to support the surface tension gradient. It was 
argued that Marangoni Flow could influence the VCHP in two ways: 

1) If the condenser end is frozen, Marangoni Flow toward that 
end would increase freezing under load. However, condenser 
freezing was argued not to be the primary mechanism for 
depriming. 

2) If the heat pipe is not frozen, Marangoni Flow will cause a 
pressure drop along fillets and natural reservoirs reducing 
their pumping capacity. However, the reduction must be 
small, otherwise anomalies due to Marangoni Flow effects 
would occur on consecutive days under similar conditions, 

not sporadically. It was argued, however, that once deprimed. 
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the open artery capacity could be measurably influenced 
by Marangoni Flow. 

In view of the fact that the control gases helium and nitrogen, par- 
ticularly helium, have solubilities in methanol decreasing with decreasing 
temperature and pressure, it was argued that this behavior could give rise 
to a potential depri.- ng mechanism which is consistent with sporadic occur- 
rences during the eclipse seasons. This mechanism is gas evolution within 
the arteries during heat pipe chilldown. 

It was postulated that when TEP ceases to operate several hours before 
an eclipse, the gas blocks the entire length of the heat pipes. During 
this off period, liquid in the arteries becomes saturated with gas surround 
ing the arteries at the gas-blocked condenser temperature and evaporator 
pressure. When the CTS enters eclipse, the heat pipes experienced rapid 
chilldown which causes the temperature and pressure to drop rendering the 
liquid in the arteries sufficiently supersaturated in dissolved gas to 
cause gas bubble nucleation within the arteries. 

In the early study, it was postulated that gas would evolve at unspe- 
cified nucleating sites leading to numerous small bubbles, some of which 
would coalesce to form fewer but larger bubbles with intrinsic longer life- 
times. When the TEP is reactivated, the heat pipes would turn on leading 
to the establishment of a gas/vapor front and a build-up of stress in the 
active portion. It was postulated then that if the local liquid stress 
exceeds a critical level for a bubble before this bubble is redissolved or 
vented, the arteries would deprime. Of the three potential mechanisms for 
artery depriming considered in this early study, gas evolution was argued 
to be the most plausible inasmuch that it ties the anomalies to eclipse 
seasons when the VCHP system experiences rapid chilldowns. Furthermore, it 
also provided an explanation for the sporadic occurrences within eclipse 
seasons since, it depends on the operating conditions before and after the 
eclipse and on bubble nucleation and aglomeration which are statistical 
phenomena . 

Although no attempts were made at that time to analytically treat the 
gas evolution scenario, a simple experiment was performed to determine 
whether gas bubbles would be nucleated in the arteries under CTS conditions 


The experimental setup consisteci of a short section of CTS-type arterial 
wick inserted into a glass tube sealed at one end and yalved at the other. 
The tube was filled with appropriate amounts of 10 percent He/90 percent 
N2 gas and methanol simulating CTS design parameters. The tube was then 
subjected to temperature reduction under hydrostatic stress in excess of 
open artery capacity. The imposed stress was equivalent to increased 
supersaturation which would enhance the probability of bubble nucleation. 
Experiments at temperature reduction rates of -1.71 C/Min and -0.9 C/Min 
yielded negative results even though the cooling rates were about twice 
those observed in orbit, and the imposed hydrostatic stress enhanced the 
probability of bubble nucleation. 

The results of this experiment lead to the conclusion that the gas 
evolution mechanism could not have caused artery depriming under conditions 
experienced by the CTS. 

Thus, th’s preliminary investigation at TRW was unsuccessful in its 
quest for a single or sequential series of mechanisms for artery depriming 
consistent with all flight data. However, depriming due to gas evolution 
within the arteries was regarded as qualitatively consistent with all the 
observed anomalies and warranted further investigation. 

1.3 PHASE I STUDIES 

In view of the fact that previous investigations of the CTS thermal 
anomalies convincingly established depriming of the heat pipes in the VCHP 
system as the cause of the anomalies but failed to identify the mechanism 
or sequential series of mechanisms that led to depriming of the arteries, 
TRW (under contract to NASA LeRC - Contract NAS3-21740) undertook a two- 
step investigation of the hydrodynamic characteristics of CTS-type heat 
pipes. 

This section of the report summarises the work performed during 
Phase I and refers to attached appendices (documents generated during 
Phase I) for further details. 

The first task performed was to identify mechanisms which under a 
wide range of feasible conditions could lead to depriming of the heat pipe 
arteries. Ten mechanisms were postulated which are briefly described and 
classified in four basic groups in Figure 1-4. Also shown are comments on 
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ARTERY DEPRIMING MECHANISMS 



[ 

\ 



- 11 - 


Figure 1-4. Observations on Postulated Artery Depritning Mechanisms 



the postulated mechanisms in light of observed characteristics of the 
anomalies and flight data. The approadh in IPhase I to inyestigate the 
potential of these mechanisms is shown in Figure 1-5. 


The basic approach has been, for a given depriming mechanism, to con- 
duct analyses in the following sequential order with the implementation of 
each succeeding task being contingent upon the results of the previous one: 

a) Analytical ";back-of-the-enyelope" calculations and/or flight 
data review. 

b) Calculations with more realistic analytical models,. 

c) Bench-scale experiments such as with glass heat pipe. 

d) Experiments with CTS SN009 heat pipe which is similar to 
heat pipe 1 on the CTS. 

e| Correlations of experimental data with thermal (SINDA) and 
heat pipe (VCHRDA) analytical models results. 

Tasks a, b, c, and d were performed by TRW, whereas task e was the result 
of a joint effort with LeRC in which the experimental data generated at 
TRW facilities were correlated with the results of steady state and tran- 
sient thermal analyses of TEP performed at LeRC. The LeRC study is briefly 
reviewed in Appendix A. 4.1. 

Mechanisms Exceeding Heat Pipe Capacity 

The arteries in a heat pipe can certainly be caused to deprime if 
the imposed heat load exceeds the designed heat transfer capacity of the 
heat pipe. Exceeding the heat pipe capacity can be the result of increased 
dissipation by TEP, excess skew in load partitioning on the VCHPS, and high 
transient heat pipe load due to condenser and/or reservior shadowing. 


Increased Dissipation (Column D of Figure 1-4) . This is clearly not 
the cause of artery depriming since: 

1) TEP power dissipation exceeding designed levels is not 
supported by f 1 i ght power data 

2) The heat loads are well below the system capacity when 
all three heat pipes are primed. In fact a single primed 
heat pipe can carry the maximum designed dissipation heat 
load. 
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Figure 1-5. ApproaGh to Investigate Potential Depriming Mechanisms 



3) On many days preceding and/or following the anomalies, higher 
heat loads were accommodated without difficulty. 

Excess Skew in Load Partitioning (Column E of Figure 1-4). The three 
heat pipes of the VCHPS on CTS are designed to turn on sequentially so as 
to balance the load between them at full-on design conditions. However, 
very low sink temperatures during equinox conditions could result in the 
load of heat pipe HPl exceeding transiently its capacity before HP2 and 
HP3 take up a significant share of the load. Failure of HPl by the above 
mechanism or by any other could then result in a rapid transfer of load to 
HP2 which, due to inertial effects associated with rapid load transfer, 
could fail leading to a rapid transfer of load to HP3, yielding, perhaps, 
a "domino effect" failure of the whole system. 

Failure of HPl due to excess skew in load partitioning is discounted 
since LeRC transient thermal results for periods preceding the four anoma- 
lies indicate the maximum load on HPl never exceeds 80 watts, a load that 
is well below its primed heat transfer capacity. 

In Appendix A.l flight data are analyzed to determine whether there 
exists a regularity in the depriming sequence associated with the anomalies, 
and the possibility of a "domino effect" failure of the system is investi- 
gated experimentally. 

The analysis of flight data showed that these exists some regularity 
in the depriming sequence as follows; HPl deprimes first, HP2 deprimes 
next followed by depriming of HP3. (On days 89 and 253 the temperature 
data suggest that heat pipe 2 or 3 was deprimed, while HPl carried the 
load.) Experimentally, attempts were made to deprime the SN009 heat pipe 
with rapid increases of load, simulating, under more stringent conditions, 
rapid sequential load transfer as heat pipes deprime. 

Tests performed for both high (-18C) and low (-96C) sink temperatures 
and at 2.5 cm evaporator tilt, indicated no measurable rate effect attribut- 
able to fluid inertia and thus argued against a "domino effect". 

High Transient Heat Pipe Load Due to Condenser and/or Reservoir 
Shadowing (Column F of Figure 1-4). Sudden shadowing cools the heat pipe 
reservoir which causes the vapor/gas front to move toward the reservoir 
end of the radiator increasing the active portion of the condenser which 
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causes the heat load on the pipe to increase. In addition, shadowing of 
the VCHPS radiator increases the load on the heat pipes due to a drop in 
the effective environment temperature. 

Although both effects will eventually lead to cooling of the heat 
source, movement of the gas/vapor front reducing the active condenser sec- 
tion, and reequilibration at the dissipation load, the transient peak may 
be sufficient to exceed the heat pipe capacity and deprime the arteries. 

An examination of flight data corresponding to periods preceding and 
including the onset of the four anomalies shows a scenario of decreasing 
temperatures of the gas reservoirs and on sections of the VCHP system 
radiator. Telemetry data for temperature sensors HPT5 and HPT6 (see Fig- 
ure 1-3) show that at these locations the temperature drop monotonically 
prior to the anomalies. 

In order to explore the potential '7f this postulated mechanism for 
artery depriming, the transient loads possible under typical orbital con- 
ditions were estimated using a simplified analytical model of the VCHP 
system. The description of the model and the basis for the transient loads 
calculations are presented in Appendices A. 2.1 and A. 2. 2, which address the 
effects of shadowing of the reservoir and condenser, respectively. The 
analytical estimates show, that under some assumed cooling rates, the 
instantaneous heat loads may be as much as 45 percent larger than steady 
state conditions would indicate. These analytical results were found sig- 
nificant enough to warrant further tests on the SN009 heat pipe. 

Appendix A. 2. 3 describes a series of tests performed on the SN009 
heat pipe during which the condenser or the gas reservoir was cooled at 
various rates while the evaporator was maintained at a constant heat load. 
In tests loads of 100 and 150 watts and cooling of the condenser at 
approximately 2.8 C/Min produced no depriming. No depriming was observed 
either in tests with 125 watts and at a condenser cooling rate of approxi- 
mately 3.9 C/Min. Thus, depriming due to rapid chilldown of the condenser 
was not demonstrated. 

During tests in which the gas reservoir was cooled repeatedly at suc- 
cessively higher rates depriming was observed. With 100 watts at the 
evaporator and the inactive condenser at -18 C, depriming occurred when the 
reservoir was cooled at the threshold rate of approximately -3.2 C/Min. 
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Although depn'ming due to rapid chi 11 down of the reservoir was 
demonstrated, the required rates were substantially higher than those 
observed during the anomalies. This fact eliminates the postulated sudden 
cooling mechanism as a potential cause of the CTS anomalies. 

Mechanisms for Depletion of Liquid Inventory 

Another possible cause of artery depriming in a heat pipe is depletion 
of liquid inventory, either depletion of total heat pipe inventory such as 
that resulting from a leak, or depletion of liquid from the active portion 
of the pipe such as could result from diffusion freezeout, suction freeze- 
out or freezing blow by. Suction freezeout is the loss of liquid from 
freezing. Freezing blowby is the loss of liquid from the thawing of an 
ice plug allowing a high pressure evaporator to blow liquid from the active 
pipe into a low pressure gas reservoir. 

Loss of inventory due to leaks will cause permanent changes in heat 
pipe performance, an event that is not supported by the observed complete 
recovery of the system following an anomaly. Diffusion freezeout is the 
transfer of inventory in the vapor phase from the active pipe into the 
frozen section where the vapor condenses and freezes. It is too slow a 
transfer process to have been effective during the time available prior to 
the anomalies. 

Suction Freezeout (Column B of Figure 1-4). Suction freezeout has 
long been recognized as an artery depriming mechanism. Depriming occurs 
due to depletion of evaporator liquid because of the local density increase 
as a freezing front moves toward the evaporator end of the pipe. It is 
necessary, however, that the condenser freeze while the heat pipe is under 
load so that the natural liquid reservoirs in the evaporator do not contain 
excess liquid. Such excess liquid would satisfy the demand of the freezing 
process without stressing the arteries to failure. Suction freezeout is a 
viable explanation for depriming on anomaly days 75, 82, and 101. But not 
on day 253, since the results of steady state and transient calculations on 
the CTS model performed at LeRC indicate no freezing on this day at the 
time of the anomaly. 

This depriming mechanism was investigated analytically in the previous 
( 2 ) 

TRW's study' ' where the mechanism is referred to as "Condenser Freezing". 
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The results of this investigation made in light of flight data and 
observation about the anomalies were summarized in Section 1.2 of this 
report. 

In Phase I of this program, suction freezeout was subjected to experi- 
mental investigation with tests on the SN009 heat pipe. Several attempts 
to deprime the arteries by freezing the condenser while the SN009 heat pipe 
was under hydrostatic and thermal load were unsuccessful. 

The experimental results support the conclusion of the previous inves- 
( 2 ) 

tigation' ' which states that although suction freezeout or condenser 
freezing can and may lead to artery depriming, it is not the primary cause 
of the anomalies. 

Freezing Blowby (Column C of Figure 1-4). The CTS heat pipes contain 
substantial amounts of excess liquid. Such excess inventory is required 
for successful priming of the arteries in earth gravity testing. In zero 
gravity a slug of excess liquid generally bridges the condenser vapor 
spaces, acting as a moving membrane between portions of the non-condensible 
control gas. Subfreezing radiator temperatures can cause the slug of 
liquid to become an immobilized plug of ice separating the active section 
of the pipe from its gas reservoir. If during a transient a pressure dif- 
ferential develops across the ice plug, with the pressure higher on the 
evaporator side, and then the plug partially thaws, this pressure difference 
can blow liquid from the evaporator to the reservoir side and deplete the 
active pipe inventory causing the arteries to deprime. 

Because freezing blowby requires freezing, this mechanism can not be 
the cause of the fourth anomaly, since freezing had not occurred immediately 
preceding the anomaly on day 253. However, freezing blowby could account 
for some of the anomalies during which freezing occurred. 

To explore experimentally the depriming potential of this mechanism, 
it was first necessary to develop means of forming an ice plug in the 
laboratory in the absence of the natural slugging of excess liquid expected 
in zero gravity. Appendix A. 3 describes a successful technique that was 
developed to form an ice plug in the SN009 heat pipe, and the procedure 
followed during the blowby tests. Three tests were performed during which 
depriming occurred each time. The results of these tests clearly established 
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that freezing blowby is a bonafide depriming mechanism and a candidate to 
explain at least some of the anomalies. 

The latter statement is supported, to a certain extent, by some results 
from the CTS model analyses which showed that ice plugs formed in all the 
heat pipes on day 82 several hours prior to the onset of the anomaly. In 
addition, the results showed that the pressure on the evaporator side 
would have been higher than on the reservoir side when the plugs thawed 
forty five minutes before the last heat pipe deprimed. The significance 
of these calculated results is that they show conditions for freezing 
blowby may have existed on at least one day of the anomalies. 

Mechanisms for Bubble Nucleation 

The presence of gas inside the arteries of the heat pipe has long been 
recognized as detrimental to the stable operation of the pipe. Because the 
solubilities of nitrogen and, particularly, helium in methanol, decrease 
with decreasing temperature and pressure, there exists the potential of 
bubble nucleation within the arteries in CTS-type heat pipes as they 
undergo rapid temperature reduction. This potential is compounded by the 
decrease in pressure that results from temperature reduction in the closed 
heat pipe environment. Liquid methanol saturated with gas can become 
supersaturated with gas when undergoing rapid temperature and pressure 
reduction, an event that enhances the potential of bubble nucleation. 

Bubbles that might nucleate inside the arteries during the transient 
cooldown of the heat pipes after the power is turned off and the spacecraft 
enters into eclipse, can coalesce to form fewer but larger bubbles. If, 
as the result of their size and the rather slow process of gas diffusion 
back into surrounding liquid, the bubbles survive until the heat pipes are 
once again activated under load, some bubbles might be convected into high- 
stress liquid regions where they can grow in size and consequently deprime 
an artery. 

The process of bubble coalescense and migration is recognized as 
statistical in nature. 

The above postulated mechanism was examined during a previous TRW 
( 2 ) 

investigation;' ' however, experiments under simulated anomaly conditions 
failed to induce bubble nucleation. 
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Despite these results, this mechanism was reexamined during the 
current program due to the fact that the postulated conditions for the 
mechanism to be activated can be supported by flight data corresponding 
to all the anomalies, and the probabilistic nature of the mechanism bears 
positive correlation to the sporadic occurrences of the anomalies. In 
addition, a thorough examination required the performance of experiments 
simulating more realistically the operating characteristics and anomaly 
conditions of the CIS heat pipes and investigation of the potential of 
other mechanisms, particularly freezing- thawing, to induce bubble nucleation. 

Bubble Nucleation Due to Temperature and Pressure Reduction (Columns G 
and H on Figure 1-4). An analysis of the potential for bubble formation 
was performed based on fundamentals described in Appendix A. 5.1. Sample 
calculations were made for two cases (1) a condenser depressurization case 
where the evaporator temperature is dropped at constant gas-blocked con- 
denser temperature and (2) a condenser chilldown case where the gas-blocked 
condenser is cooled at constant total pressure. 

The sample calculations presented in Appendix A. 5. 2 showed the possi- 
bility for bubble formation upon (1) depressurization by reduction in heat 
pipe loading and reservoir chilling and (2) condenser chilldown. Depres- 

3 

surization showed greater potential, for approximately 500,000 bubbles/cm 
were found possible, compared to only one of equal size resulting from 
condenser chilldown. The above calculations, based on two hypothetical 
cases, were followed by a number of calculations for various anomaly con- 
ditions. The predicted number of bubbles of varying critical size was 
found substantial for all anomaly conditions considered in the analysis. 

In order to verify qualitatively the above analytical results, a 
series of simple experiments were formulated which are described in Appen- 
dix A. 5.3. Tests were performed with a glass vessel half filled with 
liquid methanol containing a short section of mesh screen artery laying on 
its bottom. The vessel was instrumented to permit continuous monitoring 
of temperature and pressure. The liquid methanol was saturated with either 
helium or nitrogen at room temperature subsequent to which three series of 
tests were performed. In the first two series of tests the saturated 
methanol was subjected to various rates of temperature and/or pressure 
reduction. During the tests bubble nucleation in the bulk of the liquid 
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or on the surface of the artery mesh screen was never observed. These 
results suggested that temperature and pressure reduction are not potential 
mechanisms for bubble generations in CTS-type heat pipes. In the lest 
series of tests with the glass vessel, the saturated methanol was subjected 
to several freeze/thaw cycles. Large numbers of bubbles were observed 
streaming from the ice surface. 

Bubble Generation Due to Freezing (Column I of Figure 1-4). The 
results of the glass vessel tests were of paramount importance, for they 
identified freezing and thawing of gas-saturated methanol in the arteries 
as a potential mechanism for bubble formation in the arteries. The fact 
that freezing occurred preceding three of the four anomalies and sometime 
during the 24 hours before all of them, the involvement of ice-generated 
bubbles in the four thermal anomalies was shown to be a distinct 
possibility. 

To explore this bubble nucleation mechanism in a realistic heat pipe 
environment, a series of tests were performed with an existing glass heat 
pipe. The pipe contains a slab wick with a CTS-type artery attached to 
one side and permits observations on the behavior of the artery in an 
operational heat pipe. For the tests the heat pine was gas loaded with 
a 90 percent N 2 /IO percent He gas mixture at a pressure equivalent to that 
in the CTS heat pipes. The experiments performed on the glass heat pipe 
are documented in Appendix A. 5. 4. 

The results of bubble nucleation experiments showed that each time 
the artery in the condenser section underwent freezing and thawing bubbles 
were observed inside the arteries although the number of bubbles, their 
size and location along the artery varied from test to test. The time 
required for the bubbles to disappear by diffusion was found to vary enor- 
mously depending on their size. Small spherical bubbles dissolved within 
hours, while elongated bubbles required up to several weeks. 

Additional experiments were conducted with this heat pipe during 
v/hich the gas/vapor front was forced to move into the previously frozen 
condenser section in attempts to incorporate ice-generated bubbles into 
the active condenser and cause depriming. A number of deprimings were 
observed, but they were sporadic, probabilistic in nature. These results 
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conclusively demonstrated that; 1) Control gas is liberated from 
sarurated methanol every time the heat pipe goes through a freeze/thaw 
cycle, 2) the number, size, and durability of gas bubbles generated within 
the arteries have a statistical variation and are influenced by bubble 
coalescence, and 3) arterial bubbles can lead to depriming if they migrate 
or are convected into the active region of the pipe under normal load con- 
ditions. Thus, it is clear that freezing and thawing of the condenser can, 
but does not necessarily, lead to artery depriming, depending on subsequent 
history. 

Attempts to induce arterial failure in the SN009 heat pipe due to the 
ice-generated-bubble mechanism were partially successful. Depriming was 
not observed in several tests when the heat pipe was brought under a heat 
load immediately after the condenser had gone through a freeze/thaw cycle. 
However, depriming was discovered on two occasions after the SN009 heat 
pipe had been left unattended under hydrostatic load for several hours 
during which time the condenser froze following the termination of a normal 
freeze/ thaw cycle. 

The results of tests with the SN009 and glass heat pipes clearly indi- 
cated that in order to establish the statistics of this artery failure 
mode, repeated experiments would be required. 

SUMMARY 

In the course of Phase I of this program, numerous depriming mechanisms 
were postulated and subjected to analytical and experimental investigation. 
Several mechajiisms were convincingly rendered highly improbable causes of 
artery failure on CTS. Four mechanisms, however, were identified as poten- 
tial candidates to explain the anomalies and are listed in decreasing like- 
lihood in Table 1. 

It was argued during this Phase of the study that potential depriming 
mechanisms must be consistent with three key observations about the 
anomalies; 

1) The anomalies occur only during the radiator eclipse season 
when condenser freezing may be realized. 

2) The anomalies are sporadic. Similar conditions on successive days 
can yield an anomaly on one day and not on any of several similar 
days. 
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Table 1. Summary of Potential Depriming Mechanisms from Phase I. 


ARTERY DEPRIMING 

ANALYTICALLY 

EXPERIMENTALLY | 

RESULTS FROM: 


FLASK 

GLASS H.P 

SN-0Q9 H.P 

1. BUBBLE NUCLEATION 
DUE TO: 

a) FREEZE/THAW 


YES 

YES 

7 

b) PRESSURE REDUCTION 

YES 

NO 

NO 

— 

c) TEMP REDUCTION 

YES 

NO 

NO 

“ — 

II. FREEZING BLOW-BY 

YES 

— 

— 

YES 

III. RADIATOR/RESERVOIR 

YES 

■ 


YES 

RAPID COOLING 




IV. SUCTION FREEZEOUT 

YES 

— 

— 

NOI7J 


3) On day 253 condenser freezing did not occur immediately 
prior to the onset of the anomaly. 

In light of the above observation, the following comments can be made on 
the mechanisms listed on Table 1. 

Suction freezeout is not a statistical depriming mechanism and cannot 
account for the anomaly on day 253. Artery depriming due to suction 
freezeout was not demonstrated with the SN009 heat pipe which seems to 
indicate that none of the other anomalies could have been caused by this 
mechanism alone. This mechanism can only be considered a potential contrib 
uting factor to some of the anomalies. 

The sudden cooling mechanisms could account in principle for all the 
anomalies, however, they are not statistical in nature. Furthermore, tests 
with the SN009 heat pipe showed that the cooling rates required to fail the 
arteries far exceed those indicated by flight data. Thus, this artery fail 
ure mode can be assigned a low probability of success and only be consid- 
ered a contributing factor to some of the anomalies. 

Artery depriming due to freezing blowby was verified in repeated 
experiments with SN009 heat pipe, establishing freezing blowby as a 
bonafide potential depriming mechanism. Analyses of the CTS thermal model 
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for day 82 predicted the formation of ice plugs in the CIS heat pipes and 
calculated pressure differential accross the barrier (pressure higher on 
evaporator side) which were sustained until the ice plugs thawed. These 
predicted results tend to indicate that conditions for blowby may be 
present on some of the anomaly days. Although freezing blowby is not a 
statistical mechanism and can not be the cause of the anomaly on day 253, 
it can not be totally discounted. 

In view of the results of glass flask and glass heat pipe experiments, 
the ice-generated bubble mechanism appears to be a prime candidate for 
explaining all four CTS anomalies. It is consistent with (1) their sea- 
soned occurrence (eclipse conditions are necessary to cause condenser 
freezing), (2) their sporadic occurrence (due to the statistical nature of 
bubble population and behavior), and (3) the lack of freezing on day 253 
(bubbles were generated by a freeze/thaw process during day 252). 

1.4 OBJECTIVES OF PHASE II STUDIES 

The Phase I studies demonstrated that bubbles are generated each time 
liquid methanol containing dissolved gas undergoes freezing and thawing, 
and that these bubbles may result in artery depriming if inducted into the 
active portion of the heat pipe. This depriming mechanism was demonstrated 
on the glass heat pipe. However, depriming due to ice-generated bubbles 
was not conclusively demonstrated on the SN009 heat pipe. Because the 
SN009 heat pipe is similar in overall dimensions and performance charac- 
teristics to heat pipe number 1 onboard CTS, the need to demonstrate 
depriming on the SN009 due to bubbles became apparent in order to demon- 
strate beyond a reasonable doubt that the postulated bubble mechanism is 
the cause of the thermal anomalies. 

The Phase I studies also warranted further study of the behavior of 
bubbles inside arteries, particularly their lifetimes, owing to the fact 
that the probability of arterial failure due to bubbles induction into 
•^igh stress regions is positively correlated to longer bubbles lifetimes. 

As a result a second phase of this program was continued. The second 
phase was necessary to bring the CTS-anomaly studies to the satisfactory 
conclusion envisioned at the outset. 
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The objectives of Phase II were: 


1.. To analyze the idiffusion iprocess for non-condensiibTe control 
gases in arteries and to determine lifetimes of bubbTes of 
various sizes at yarious temperature levels and rates of 
change of 'temperature level and for various selected control 
gases,. 

2. To run continuous cyclic freeze/thaw tests on the SN009 heat 
pipe for fiye consecutiye days in an effort to produce artery 
depriming by freeze-thaw ibubble formation. The successful 
triggering of a failure would establish that freeze-thaw bubble 
formation i,^ the logical cause of the CTS thermal anomalies. 

3. To propose concepts to avert artery depriming. 
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ZjQ ftNALYSIS OF BUBBLE LIFETIMES IN HEAT PIPE ARTERIES 


2.1 THE SPHERICAL BUBBLE MODEL 

Consistent with Phase II objectives 1 and 3, It Is desired to calculate 
how fast vapor bubbles in heat pipe arteries might grow or collapse and 
how the !ki;nd of gas used affects the growth/collapse rates. The first 
imodel selected for such calculations is the one-dimensional spherical 
annulus, the domain R^ < r < Rg as shown in Figure 2-1. 



Figure 2-1. Spherical Bubble Model 

The vapor-and-gas-filled bubble (r < Rj) has uniform composition, 
hecause gas-phase diffusivities exceed liquid phase diffusivities by 
several orders of magnitude. The external gas r > R 2 has known pressure 
and composition. The interface at R 2 is considered to be a screen like 
the wall of an artery so that the pressure difference across the interface 
is uncoupled from that otherwise dictated by surface tension and R 2 . At 
r=R| no screen exists, so the vapor-to-liquid pressure difference is 2a/R| 
where a is surface tension. 
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The foil owing heat pipe scenario is postulated: The heat pipe 
pictured in Figure 2-2 operates with an active section at 300K and a gas- 
blocked condenser at 250K long enough for the artery liquid to saturate 
with the gas. Then as pictured in Figure 2-3 the condenser is briefly 
frozen and then thawed to 18DK. Upon thawing of the gas-blocked condenser, 
there are, as has been observed experinientally, large numbers of minute 
bubbles released from the ice, but most diffuse back into solution. A few 
are postulated to condense to form the bubble of prescribed size R 2 q. 

The noncondensible gas compositions in the bubble and in the liquid were 
taken to be the same, the same as that in the saturated liquid at 250 K. 

The condenser then warms up at a prescrsbed rate as pictured in Figure 2-3, 
affecting the vapor composition in the external gas. Because thermal 
diffusion proceeds at a much higher rate than mass diffusion (Prandtl num- 
ber is much smaller than Schmidt number), the artery liquid is assumed to 
be in thermal equilibrium with the condenser temperature, and the vapor 
composition of the internal gas is also affected. 
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Figure 2-2. Axial Temperature Profile 
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Figure 2-3. Postulated Gas-Blocked Condenser Temperature History 


2.2 GOVERNING EQUATIONS FOR THE SPHERICAL BUBBLE 

Within the liquid, conservation of mass species i takes the form. 



+ V 


9X. 


9r 



(2-1) 


where i = 2 to n (i=l is reserved for the liquid). The quantity is the 
mole fraction of species i, and r is the radius. The assumption has been 
made that the solution is dilute and isothermal, so the molar concentration 
c is constant at that of the liquid, and the diffusivity for species i 
in the liquid is unaffected by the other dilute species present. The 
radial velocity v is caused entirely by growth or collapse of the central 
bubble, because the diffusion mass transfer rates are so low, conservation 
of mass dictates 


V = 



dRl 

dt 


( 2 - 2 ) 
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Within the central bubble the gas is assumed to be uniform in 
composition and pressure, because the gas-phase diffusivities are high 
compared to the liquid-phase values, and inertia forces are negligible for 
the small changes in bubble growth or collapse rate. Accordingly the gas- 
phase species equations equivalent to Eq. (2-1) are not used explicitly. 
Rather, they are replaced by the simple expression 


dw- 

■Jt 


i = -4irR? N 


1 ”i,l 


(2-3) 


where w. is the molar content of species i within the bubble and N.- , is 

1 I j 1 

the molar flux for species i at r=Rj, given by Pick's Law on the liquid 
side of the bubble meniscus as 


3x. 

N. , = c D. - 7 -^ 
1,1 1 8r 


r=R, 


i=2,n 


(2-4) 


O 

where c is liquid molar concentration (g-moles/cm ). 
The molar content of the bubble is given by 


w. Rj^ P./R i=2,n (2-5) 

where is the condenser temperature, and R is the universal gas constant. 
Similarly the conservation of momentum equation is replaced by a quasi - 
hydrostatic balance equation 


E P, = P„q + 2 c/Rj (2-6) 


The partial pressures within the bubble are given by Raoul t's and 
Henry's Laws. Since is nearly unity, Raoult's Law gives simply 






(2-7) 


while Henry's Law is 


Outside the liquid at r=R2 the total pressure of the vapor and gas is 
Pv(Tgv) where T^^ is the evaporator temperature. The gas pressure in the 
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gas-b1ocked condenser Is 

Pg = P, (Te,) - Pv (Tc) (2-9) 

The gas composition is specified by Y^. , i=2 to n, where Y^. is the mole 
fraction of species i in the noncondensible. For example, if 2 is helium 
and 3 is nitrogen, Y 2 might be 0.10 and Y^ 0.90 for a 10 percent helium, 

90 percent nitrogen gas mixture. The saturation values of x- at r = R 2 are 

*i = ''i Pg/6,- (Tc) (2-10) 

The initial conditions at t=0 are specified by a uniform set of x-, 
i=2 to n, for the initial liquid composition, a set of Y^. in the gas bubble 
and the initial radius of the bubble. From these there can be derived 

the set of w., i=2 to n. First, from Eqs. (2-6) and (2-7) 

'■g ' ’’liq (T,)/Ri - P„ (T^) (2-11) 


Then from the gas law 


= T TfR? Pn/RT, 
g 3 1 g c 


where 


Finally 


n 

Wq = E Wi 

9 i=2 ^ 


«i = h ^ 


2.3 TRANSFORMATION AND NUMERICAL SOLUTION 


(2-12a) 


(2-12b) 


(2-13) 


It is convenient to anchor the spatial coordinate to the bubble inter- 
face and make the coordinate dimensionless with respect to R 2 - R^. 


r - R^(t) 

" R2(t) - R^(ty 


(2-14) 
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where R 2 (t) is given by 


4 

3 


TT 





3 

1.0 


(2-15) 


Velocity v given by Equation (2-2) becomes 


V = 



(2-16) 


Partial derivatives with respect to r at constant t (and thus at constant 
Rj and AR = R 2 - Rj^) are given by 


9x^. 3x^- 


9r 


1 

9ti AR 


(2-17) 


The partial derivative with respect to t at constant r must be transformed 
to one with respect to t at constant n 


9t 


9n 

9t 


9X^. 


9X^. 


9n 

t 'r 


1 ' (’^”*^1) ’ 

- Ri - AR 

AR 1 ^^2 


Hence Equation (2-1) becomes 


(2-18) 

(2-19) 


9t 


9X 


" AR '^1 ■ AR ^7 \R^+tiAr) AR W 


2 

= D + — 2 

^ \aR^ 9ti^ '^1 I 


Collecting terms gives 


9x. 
1 ^ 

9t 


2 ' 
V «i 


^(R^+tiAR) AR 


(l+n)Rj 

a1 


2D. 

1 


9X. 
1 


Di 9 X. 


(Rl+nAR)AR/ 9n “ 3 ^^ 


= 0 


( 2 - 20 ) 


Since large values of 9x^./9ti are expected near n = 0, computational 
efficiency is improved by transforming the n coordinate to z so as to 
expand the region near n = 0 relative to the region near n=l. A negative 
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value of the parameter y in the following coordinate transformation 
achieves the desired result. 


^Yn 


z = 


1 


After transformation Equation (2-20) is in tie form of 





= 0 


( 2 - 21 ) 


( 2 - 22 ) 


where F and G are functions of z, Rj, Rj, D^, and y (recall that R 2 and 
hence AR is given by Equation (2-15)). 

Equation (2-22) is solved numerically in finite difference form. Let 
subscript j denote z location and superscript 0 and 00 denote values at 
t-At and t-2At respectively. The species subscript i is dropped to avoid 
confusion. The temporal derivative is approximated as 


^ = ^[a 
at At '■®l 


X. + a, x.° + a~ x.°°] 
J «j 


(2-23) 


where for constant time increment At, a^ = 3/2, a 2 = -2, and a^ = 1/2. 
The spatial derivatives are 


az 2 az 


(2-24) 



’‘.HI * ’’.i-i - 


Equation (2-22) takes the form 


"j 




'j+1 




'j-1 


+ C., 


j = 2,N 


(2-25) 


(2-26) 


where z« = 0 and z« = 1. The coefficients A., B., and C. may be calculated 

^ ^ n 00 

at any time step in terms of the previous known sets xV and x. and the 

• J J 

known values of R, , R, , D., y» and z.. The new set of x. values is computed 

XXI J J 

by means of Gaussian elimination in which 
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(2-27) 


ic ic 

X. = A. X.., + B. 

J 0 J+1 J 


Equations (2-26) and (2-27) combine to permit the forward (j=2, 3, ...) 

ic ic 

calculation of A. and B. . 

J J 

A 2 * = 0 , Bg* = Bg x^ + C 2 (2-28a,b) 

Aj* = Aj/(1 - Bj (2-28c) 

Bj* ' (“j 

Then Equation (2-27) allows the backward (j=N-l, N-2, ...) calculation of 
Xj starting with the prescribed boundary value x^^. 

Appendix B.l contains the listing of the program used to calculate 
x(z.,t) and R, (t). 

»J 

2.4 THE ELONGATED BUBBLE MODEL 

When a bubble grows so that its radius (Rj^) reaches the artery radius, 
further growth occurs through elongation of the bubble within the artery. 

In the absence of a priming foil, the artery wall prevents meniscus coales- 
cence, and a sheath of liquid is retained about the elongated bubble. 

Mass transfer of species i (i=2 to n) can occur through the sheath and into 
or out of the end-cap liquid. Equations (2-5) and (2-12) become 

Wi = (| TT R^ + ttR^L) P./RT^ (2-29) 

“g = <1 " "a ^ .r2l) P^/RT^ (a-30) 

while Equations (2-6) through (2-11) and (2-13) remain unchanged. 

Equation (2-3) and (2-4) are well approximated by 


dw • 

^ = CD. [(2irRgL s/tAR) + (F/R^)] Ax. (2-31) 

where Ax^ is the difference in mole fraction between that of the inside 
liquid at the bubble-liquid interface and that of the outside liquid at 
the condenser-gas-liquid interface. The quantity e is the artery screen 
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void fraction, t is its tortuosity, and F is the conduction shape factor 
for a hemispherial bubble of radius R, in a semi -infinite cylinder of 

a 

radius R, + aR. Fcr the CIS heat pipe arteries R. = 0.0127 cm, e = 0.37, 
R./AR = 8, and F = 8. 

a 

Equation (2-31) was integrated numerically with Equation (2-30) used 
to find L in a Runge-kutta type algorithm. 

2.5 PROPERTIES USED IN CALCULATIONS 

Figures 2-4 (Ref. 3) shows the solubility data used for the calcula- 
tions. Shown versus temperature is the mole fraction in the liquid for a 
partial pressure of 1 atm, that is the reciprocal of Henry's constant in 
atm'^. The variation with temperature will be shown to be quite signifi- 
cant, For now, merely note that methane and argon become less soluble in 
methanol with increasing temperature (as does air in water) while helium 
becomes less soluble with decreasing temperature. The solubility of nitro- 
gen in methanol is nearly independent of temperature. Helium is seen to 
be only sparingly soluble while methane is quite soluble in methanol. 

Figure 2-5 (Ref. 4) shows the diffusivity of argon and helium in 
liquid methanol versus reciprocal temperature. The theory shown in the 
figure agreed poorly with experimental data, and it was decided to fit the 
experimental data with an Arrhenius relation. Since data were lacking and 
theory seemed dubious, it was decided to use the experimental helium fit 
for helium diffusivity and the experimental argon fit for the diffusivity 
of argon, nitrogen, and methane. Thus the calculated results for the 
latter two gases should be regarded as only qualitatively correct. When 
reliable data for these latter two gases become available, a minor time- 
scale expansion or contraction will be necessary to adjust for the approxi- 
mate values of diffusivity used. 

2.6 RESULTS FOR THE SPHERICAL BUBBLE 

Figures 2-6 through 2-9 show the calculated results for the spherical 
bubbles. Each of the four figures is for a different gas: helium, argon, 
methane, and 10 percent He 90 percent N 2 » respectively. For a given con- 
denser warmup rate and initial bubble size, helium bubbles are seen to 
persist for very much longer times than argon bubbles. Methane bubbles 
disappear much faster than even argon bubbles. The explanation is simply 
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Figure 2-6. Spherical Bubble Radius History, Helium Gas 




Figure 2-7. Spherical Bubble Radius History, Argon Gas 






Figure 2-8. Spherical Bubble Radius Hi story * Methane Gas 



Figure 2 - 9 . Spherical Bubbl 


that, despite its higher diffusivity, the solubility of helium in methanol 
is low; hence the values of cD^. Ax^. are low for helium. Because of the 
helium component, the 90%N2 - 10%He noncondensible gas mixture displays 
long bubble lifetimes. 

For a given gas and condenser warmup rate, small bubbles disappear 

much faster than large bubbles. In the small bubble the gas pressure is 

larger due to surface tension, so the values of x^. are larger. The mass 

transfer coefficient area product decreases with size as R-i, but the 

3 ^ 

volume of gas remaining decreases as . 

With helium gas in a bubble of a certain initial size, increasing the 
condenser warmup rate increases the solubility and diffusivity of the 
helium and hastens the reabsorption of the bubble, despite a decreased 
surface tension, and an increased pressure of methanol vapor which tends 
to swell the bubble. 

With argon the increase in diffusivity with increasing temperature 
causes the same type of behavior as in the case of helium but to a lesser 
degree. With methane, however, the increase in diffusivity offsets the 
decrease in solubility and lessened surface tension, so that the time to 
reabsorb small bubbles is nearly independent of condenser warmup rate. A 
complex behavior is seen for the large bubble and high condenser warmup 
rate. In that case, the fall in solubility of methane and the increase in 
methanol vapor pressure combine to cause the methane gas and methanol vapor 
bubble to swell before ultimately collapsing. 

Table 2-1 summarizes ths bubble lifetime results. Note the much 
longer times necessary for helium-containing bubbles to be reabsorbed. 

2.7 RESULTS FOR THE ELONGATED BUBBLE 

Figures 2-10 through 2-13 show calculated results for the elongated 
bubbles. Again each figure is for a different gas or gas mixture. Again 
helium-containing bubbles are seen to persist to a much longer time, while 
methane bubbles are most readily reabsorbed. Again an increase in conden- 
ser warmup rate strongly hastens the reabsorption of helium bubbles and 
weakly hastens the Veabsorption of argon and methane bubbles. Bubble 
swelling, elongation this case, is seen in all cases, because reabsorp- 
tion is initially slower than for small spherical bubbles. A final 
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Table 2-1. Calculated Spherical Bubbles Lifetimes 




Condenser Warm-up Rate (K/sec) 

Gas/Gas 

Initial 

Radius 

(cm) 

0.0 

0.01 

0.05 

Mixture 

Time, in seconds, at which 
bubble radius reaches 1/10 
of initial radius. 

Methane 

CH 4 

0.01 

15 

15 

15 

0.02 

68 

68 

68 

0.04 

3775 

4362 

3712 

Helium 

0.01 

10717 

4055 

1658 

0.02 

64411 

9998 

6782 

He 

0.04 

371408 

38651 

35409 

Argon 

Ar 

0.01 

208 

204 

194 

0.02 

3614 

2790 

1884 

0.04 

27890 

10740 

8460 

CTS Mixture 
90% N 2 / 10 % He 

0.01 

1514 

1216 

782 

0.02 

10777 

4950 

2823 

0.04 

67520 

16038 

14150 
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Figure 2-10. Elongated Bubble Length History, Helium Gas 



TIME (SEC) 10 

Figure 2-11. Elongated Bubble Length History, Argon Gas 





observation is that the helium in the helium-nitrogen gas mixture becomes 
the residual gas in the bubble at long times and low condenser warmup rates 
and requires a long time to diffuse out of the bubble. 

labile 2-2 summarizes the calculations of the elongated bubble lifetimes. 


Table 2-2. Calculated Elongated Bubbles Lifetimes 




Condenser Warm-up Rate (K/sec) 

Gas/Gas 

Initial 

Length 

CM 

0.0 

0.01 

0.05 

Mixture 

Time, in seconds, at which bubble 
length reaches 1/10 of initial 
length. 

Methane 

CH4 

2,5 

110,070 

42,655 

40,55B 

Helium 

He 

2.5 

i 

5,069,500 

436,000 

432,300 

Argon 

Ar 

1 

2.5 

406,500 

101,600 

99,000 

CTS Mixture 
90% N2/10% He 

2.5 

998,800 

182,560 

186,440 


2.8 SIGNIFICANCE OF THE ANALYTICAL/NUMERICAL RESULTS 

The CTS heat pipes contained 90%N2 - 10%He control gas. The helium 
was added to facilitate leak detection. In retrospect, that decision 
appears to have been unfortunate. For should gas bubbles be created by 
freezing in the gas-blocked condenser and should they coalesce into bubbles 
approaching. 0.8 mm in diameter, their lifetimes are seen to be on the 
order of 20 hours when the condenser warmup rate is slow. Larger sausage 
bubbles are seen to persist hundreds of hours, that is, some ten days or 
more. Such bubbles, if present through cyclic freezing and thawing may be 
enlarged through a fractionation process reminiscent of the making of New 
England applejack. With each freezing episode dissolved gas is forced out 
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of the ice into the bubble. A large bubble can persist many days even 
after freezing no longer occurs provided it is in the gas-blocked region. 
When the heat pipe is powered up, and the gas front approaches the bubble, 
Marangoni Flow brings the bubble into the active condenser, and condensate 
flow sweeps it up the pipe toward the evaporator to a location where the 
local vapor- liquid pressure difference exceeds the capillary pressure of 
the artery, and the artery deprimes. 

The results suggest that methane would make a good control gas, 
because methane bubbles up to 0.8 mm in diameter are reabsorbed in approxi 
mately one hour. Thus there is little likelihood that a 24 hour periodic 
freezing and thawing could create a large bubble, and the small bubbles 
released from thawing ice would be reabsorbed before the approach of the 
gas front. 
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3.0 CYCLIC FREEZE/THAW TESTS ON SN009 HEAT PIPE 

3.1 TEST OBJECTIVES 

Phase II objective 2 was to simulate an anomaly with the SN009 heat 
pipe by cyclic freezing and thawing of the condenser during a 5-day extended 
continuous test period. The goals of the tests were (a) to demonstrate the 
generation and collection of bubbles in the arteries and their migration 
or induction into the active section of the heat pipe where their presence 
could result in depriming of at least a single artery, and (b) to establish 
by test data analysis that freezing blowby or suction freezeout rather than 
bubble formation was not the cause of the observed heat pipe failure. 

3.2 APPARATUS AND TEST PROCEDURE 

The test assembly described in Figure 3-1 is similar to the one used 
for tests on the SN009 heat pipe during Phase I of this program. As shown, 
the heat load is applied on the heat pipe over a 0.305 meter long section 
by means of an attached 4.0 kg heated aluminum block. 

The heat pipe condenser, 0.91 meter long, is mounted on an aluminum 
plate whose area and thermal mass corresponds approximately to the effective 
portion of the CTS radiator. This plate is coupled to the sink through 0.32 
centimeters of cork over a 0.1045 square meter area. Tape heaters attached 
along the plate allow changing the condenser temperatures even with a con- 
stant temperature sink. 

A cooling coil with tape heaters is attached to the gas reservoir for 
purposes of controlling its temperature independently from the rest of the 
heat pipe assembly. 

The heat pipe test assembly is instrumented with seventeen temperature 
sensors, the location of which are shown in Figure 3-2. The sensors are 
connected to a 24-channel strip chart recorder which allows monitoring these 
temperatures at 1.5 minute intervals. 

The test procedure in Table 3-1 outlines the basic steps for conducting 
the cyclic freeze/thaw tests on the SN009 heat pipe. Steps 1 through 5 
represent the original procedure for repriming the heat pipe and adjusting 
the operating parameters in preparation for the commencement of a new cycle. 
As indicated in the procedure, the selected evaporator elevation and heat 
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Figure 3-1. Engineering Sketch of SN009 Heat Pipe Test Assembly 











Table 3-1. SN009 Heat Pipe Freeze/Thaw Cyclic Test Procedure Outline 


Using the test assembly described in SK80010 and SK80011, perform the 
cyclic tests for a period of 120 continuous hours in accordance with the 
following procedure: 

1. Level heat pipe to within 1.0 mm. 

2. Isothermalize heat pipe to within 1°C and wait 30 minutes. 

3. Elevate evaporator end 1.2 cm above condenser. 

4. Apply 130 watts to evaporator heater and verify artery priming. 

If arteries are primed, proceed with Step 5, otherwise, repeat 
Steps 1 through 4. 

5. Reduce sink plate, inactive condenser and gas reservoir tempera- 
tures to the following levels: 

a) Sink plate: lower than -145°C (-229°F) 

b) Inactive condenser: -40°C ±5°C 

c) Gas reservoir: -95°C ±5°C (-139°F) 

Reservoir temperature reduction rates should not exceed 
-2°C/min (-3.6°F/min). 

6. Simultaneously, reduce temperatures in inactive condenser section 
below -102°C (-152°F) and increase gas reservoir temperature to 
-40°C ±5°C over a period of approximately 25 minutes. Maintain 
subfreezing condenser temperatures for at least 12 minutes. If 
heat pipe fails, repeat Steps 1 through 6. 

7. Increase inactive condenser temperatures to -40°C ±5°C over a 
period of approximately 20 minutes and simultaneously start 
reducing gas reservoir .temperature to -95°C over a period of 
180 to 200 minutes. At no time should the reservoir cooling 
rate exceed -2°C/min. If heat pipe fails, repeat Steps 1 
through 7. 

8. Repeat Steps 6 through 8. 
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load oOv the heat pipe for the tests are 1.25 centimeters and 130 watts, 
respectively. 

The tilt on the pipe is meant to minimize the contribution of excess 
fluid inventory to the pumping capacity of the pipe and to compensate to a 
certain extent for that component of the buoyancy force vector tending to 
keep the bubbles against the upper walls of the arteries. 

Although heat loads over 90 watts are known to be in excess of the 
heat pipe one-artery capacity, the higher heat load of 130 watts is chosen 
in order to enhance the migration of bubbles toward the evaporator by vir- 
tue of higher fluid flow rates inside the arteries. 

Steps 6 through 8 of an idealized freezing/thawing test cycle are 
shown graphically in Figure 3-3 where segments (B-C), (B-D), and (B-E) 
represent different condenser warmup rates. These steps are described in 
conjunction with the figure in what follows. 



Figure 3-3. Idealized Freezing/Thawing Cycle "A" 
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By turning off the heaters on the condenser plate, the inactive 
condenser section is allowed to cool from -40°C (Point A) to below the 
freezing point of methanol (B) at a rate of approximately -4°C/min. 
Simultaneously, heat is applied to the gas reservoir to raise its tempera- 
ture from about -95°C to -40°C in approximately 25 minutes. During this 
period of the test cycle, the advancement of the gas front toward the 
evaporator resulting from reduced condenser temperatures is further enhanced 
by the increasing temperature of the gas reservoir, the net effect of the 
above being to inactivate temporarily a portion of the condenser to allow 
the liquid in the arteries to freeze. Subfreezing temperatures in the 
inactive condenser are maintained for ten minutes to ensure complete freez- 
ing inside the arteries. 

Failure of the heat pipe during condenser freezing could be attributed 
to be the result of either suction freezeout or the bubbles mechanism. If 
this event is repeatable in every cycle, then artery depriming will be, 
with a high degree of certainty, due to suction freezeout. On the other 
hand, sporadic failure of the heat pipe during condenser freezing will be 
consistent with the recognized statistical nature of the bubbles mechanism. 
If no evidence of heat pipe burnout is observed through the end of the 
freezing process, the test cycle continues with Step 7, otherwise. Steps 1 
through 6 are repeated. 

During Step 7, (e.g., B-C) power is applied to heaters on the conden- 
ser plate to thaw the inactive condenser and raise its temperature in 
approximately 25 minutes to -40°C where it is held within +5°C for the 
remainder of the 240-minute test cycle. About the time the inactive con- 
denser starts to thaw (B+), cooling of the gas reservoir is initiated and 
allowed to continue through the end of the cycle (F) at an average rate of 
0.3°C/min. Reservoir cooling rates are not permitted to exceed 2°C/min to 
preclude the sudden cooling mechanism from taking place. 

During this period of the test cycle, (e.g., B-C) rising condenser 
temperatures will cause the gas/vapor front to advance toward the reservoir 
end of the condenser, thawing the previously inactive section of the pipe 
and presumably engulfing into the active section bubbles that were gener- 
ated during the freeze/thaw process. The simultaneous and continuous 
reduction of the gas reservoir temperature will have t';o significant effects 
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on the bubbles mechanism. First, it will cause a continuous movement of 
tie gas front toward the reservoir during which additional bubbles will be 
induced into the active portion of the heat pipe where their migration 
toward the evaporator is enhanced. Second, reservoir cooling will tend to 
diminish the pressurization of the heat pipe resulting from increasing 
condenser temperatures and will cause subsequently, after the condenser 
reaches a steady state, (e.g., C) a continuous decrease of the total pres- 
sure in the heat pipe. 

Exploratory analyses using the spherical bubble model indicate that 
pressure reduction in the heat pipe environment can significantly increase 
bubbles lifetimes. Longer bubble lifetimes are certainly a factor enhanc- 
ing the potential of the bubble mechanism to deprime the arteries. 

A heat pipe failure observed about the time the inactive condenser 
undergoes thawing (B+) could be attributed to either the postulated freez- 
ing blowby or the bubbles mechanism. The fact that artery depriming due 
to freezing blowby is not statistical in nature requires in order to con- 
sider the freezing blowby mechanism the culprit, that failures during the 
thawing process (B+) occur at every cycle. Sporadic heat pipe failures 
during the thawing process or at any other time during the test cycles will 
be considered to be the result of bubble migration and growth in the 
arteries. If no evidence of heat pipe burnout is observed by the end of 
the 240-minute test cycle, (F) Steps 6 and 7 are repeated for the remainder 
of the 120-hour continuous test period. 

3.3 TEST RESULTS 

The cyclic tests on the SN009 heat pipe were performed from April 14 
to April 19, 1980. A total of twenty seven partial and complete cycles 
were achieved during which seventeen burnouts or heat pipe failures were 
observed. A summary of the test results is presented in Figure 3-4. It 
shows the temperatures at nine different locations along the heat pipe at 
the end of the cooling period (Condition 'A') and those observed near the 
time of burnout or at the end of a cycle (Condition 'B'). 

The tests commenced at 0600 hours when the heat pipe was tilted 
1.25 cm and 130 watts were applied to the evaporator block. After having 
esLablished that both arteries were primed, cooling of the condenser was 
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initiated by turning the radiator heaters off. This process continued 
until subfreezing temperatures in the inactive condenser section, as indi- 
cated by temperature sensors 10, 11 and 12, were realized. These sensors 
showed minimum temperatures of -120C, -132C and -114C respectively prior 
to the initiation of the condenser warm-up period. Approximately eleven 
minutes after the radiator heaters were turned on, the heat pipe failed. 
About the time of the failure, the condenser warm-up rate (CWR) was about 
9F/min and the minimum temperature in the condenser section, shown by sen- 
sor 11, was -79C which is 19C above the freezing point of methanol, thus 
well past point in Figure 3-3. 

Subsequent to burnout, the heat pipe was leveled and allowed to reach 
close to isothermal conditions over a period of two hours. The pipe was 
then tilted 1.25 cm and 130 watts were applied to the evaporator, but the 
heat pipe failed to sustain the load. The above priming procedure was 
repeated with similar results. 

The inability of the heat pipe to hold the load after two priming 
attempts was postulated to be the result of the presence of residual 
bubbles in the arteries. It was decided then to modify the procedure to 
prime the pipe, by requiring that the evaporator end be raised 61 cm for 
one minute before leveling the pipe. The high tilt would presumably eject 
all the liquid and residual bubbles from the arteries. After the pipe was 
leveled for thirty minutes, 130 watts were applied to the evaporator raised 
1.25 cm. This time priming of the arteries was verified and the second 
cycle was started. The latter procedure for priming the heat pipe was 
followed through the rest of the test period with successful results. 

Cycle number 2 was the first complete normal cycle and was followed 
by three burnouts observed during the condenser warm-up periods of cycles 3, 
4 and 5, at which times the CWR was, approximately 5 C/rain (9F/min). 

The partial history of selected temperatures along the heat pipe dur- 
ing normal cycle 2 is shown in Figure 3-5. It can be seen that prior to 
condenser warmup the ice front almost reaches the locations of sensor 9, 
i.e., approximately 18 inches of the condenser end were frozen. During 
the warmup period of the cycle the figure shows the gas front moved close 
to sensor 10 which indicates that a portion of the previously frozen con- 
denser became part of the active condenser section. 
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TIME (MINUTES) 

Figure 3-5. Temperature Histories During SN009 Heat Pipe Cyclic Test No. 


The large proportion of failures, four out of five, was attributed to 
rapid thawing and warming of the inactive condenser, events which were 
originally postulated to result in the formation of bubbles inside the 
arteries and in their rapid induction into the active portion of the con- 
denser before the bubbles could vent by diffusion. Once in the active 
portion of the pipe, these bubbles would grow and deprime the artery. 

A typical temperature history of the four cycles during which deprim- 
ing was observed is that of cycle 4 shown in Figure 3-6. As shown, the 
ice front was located somewhere between sensors 9 and 10 at the end of the 
cooling period of the test cycle. During the warmup period depriming is 
observed apparently at the time the gas front reached the previous location 
of the ice front. The fact that depriming occurred approximately 9 minutes 
after the inactive condenser completely thawed, the possible involvement of 
freezing blowby in these heat pipe failures can be discounted, therefore 
it can be argued that depriming can only be the result of bubbles. 

In order to determine the effect of condenser warm-up rates on the 
frequency of heat pipe failures, the CWR during the next six cycles was 
reduced an average of fifty percent. The results were that the heat pipe 
performed normally during two consecutive complete cycles, cycles 6 and 7, 
at the reduced CWR of about 2.2C/min. The following four cycles, however, 
resulted in burnouts during condenser warmup at rates varying from 2.2C/min 
to 3.6C/min. The variation in these CWRs was not intentional but rather 
the result of the limited heate>^ control capability of the test assembly. 

The typical partial temperature history of normal cycles 6 and 7 is 
that of cycle 6 shown in Figure 3-7. The temperature history of cycle 11 
shown in Figure 3-8 is similar to that of cycles 8 and 10. As in most 
previous cycles, the ice front can be seen located between sensor 9 and 10 
at the end of the cooling period. Figure 3-8 shows however, that depriming 
occurred under significantly different condenser conditions in most cycles 
at the lower condenser warmup rate: 1) the inactive condenser was still 

partially frozen (last 12 inches) and 2) temperatures of sensors 8 and 9 
indicate the gas/vapor front never reached the previously frozen region. 

Because the condenser was not completing the thawing process at the 
time of burnout, freezing blowby can be discounted as the cause of pipe 
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failure which, then can be attributed only to be the result of bubbles. 
However, the fact that portions of the previously frozen condenser section 
were never incorporated into the active heat pipe section, a condition that 
was postulated to be necessary to induce bubbles into the evaporator, 
seems to indicate that other mechanisms such as Marangoni Flow were involved 
to move bubbles into the active portion of the condenser where further 
growth or movement would be expected. 

It was postulated that depriming due to bubbles released in the inac- 
tive section occurred as follows. During the cooling period the freezing 
front moves toward the evaporator expelling bubbles which are then locked 
inside the ice and building stress on the arteries due to suction freeze- 
out. This process continues until the end of the cooling period when the 
imposed stress on the arteries due to freezeout reaches a maximum level. 

This increase in stress is accompanied by a reduction of the hydrodynamic 
stress which results from shrinkage of the effective length due to advance- 
ment of the gas/vapor front toward the evaporator. 

This imposed stress on the arteries due to suction freezeout is either 
overcompensated by the reduced hydrodynamic stress or the unused pumping 
capacity of the artery, for no heat pipe failures were observed during the 
cooling periods of the test cycles. During the warmup period, the advance- 
ment of the gas/vapor front toward the reservoir end results in increased 
hydrodynamic stress as the effective length expands, and in the propaga- 
tion of a conduction heat wave which causes the ice front to recede and 
release liquid and bubbles upon thawing. Although the release of the 
liquid reduces the stress that came from suction freezeout, there remains 
sufficient back stress on the arteries in the gas-blocked region to induce 
growth of the released bubbles. These bubbles agglomerate and then reach 
the active section due to continuing growth and/or Marangoni Flow. 

Whereas depriming during cycles 8, 10, and 11 occurred 10 to 12 min- 
utes into the warmup period before the gas/vapor front reached previously 
frozen sections, depriming in cycle 9, as shown in Figure 3-9, occurred 
approximately 45 minutes into the warmup period, 25 minutes after the con- 
denser completely thawed. An interpretation of the data seems to indicate 
that bubbles first released about the previous location of the ice front 
somewhere between sensors 9 and 10 were subsequently incorporated into the 
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Figure 3-9. Temperature Histories During SN009 Heat Pipe Cyclic Test No. 




active portion and caused depriming when the gas/vapor front reached that 
location 45 minutes later. 

During test, cycles 12 through 18, the heat pipe inactive condenser 
was warmed up at even lower rates ranging from 0.4C/min to 1.7C/min. 

Burnouts were observed in all these cycles, except in cycle 18 which had 
the lowest CWR of 0.4C/min. 

The temperature history of cycle 13 is shown in Figure 3-10. It can 
be seen that depriming occurred before the gas/vapor reached previously 
frozen sections and while the last 12 inches of the condenser were below 
the freezing point. Similar observations can be made of deprimings 
observed durinvg cycles 12, 14, 16 and 17. 

Depriming during cycle 15, as shown in Figure 3-11, occurred under 
different conditions. It can be seen that when the gas/front apparently 
reached previously frozen areas a burnout was triggered. At this time the 
condenser had been completely thawed for approximately 40 minutes. 

Normal freezing cycle 18, was followed by three cycles during which 
the inactive condenser was rapidly cooled to temperatures just above the 
freezing point and then warmed at rates as high as 6.7C/min. The idealized 
nonfreezing cycle is described graphically in Figure 3-12. These nonfreez- 
ing test cycles were performed to support *he hypothesis that the frequency 
of failures due to bubbles could be greatly diminished if freezing and 
thawing did not occur in every cycle. The fact that the heat pipe evidenced 
normal operation during these cycles tends to support the above. The typi- 
cal temperature history of a nonfreezing test cycle is that for cycle 19 
shown in Figure 3-13. 

During test cycles 22 through 24, the inactive condenser was rapidly 
cooled below the freezing point with the heat pipe operating under 130 watts. 
Subfreezing condenser temperatures were maintained for at least ten minutes 
subsequent to which the evaporator heaters were turned off and the inactive 
condenser was rapidly warmed. When the inactive condenser temperatures 
reached about -40C the evaporator heaters were turned on again. The 
idealized cycle is described graphically in Figure 3-14. A burnout was 
observed during the second cycle, cycle 23, 15 minutes after the heat load 
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Figure 3-11. Temperature History During SN009 Heat Pipe Cyclic Test No. 15 
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Figure 3-13. Temperature History During SN009 Heat Pipe Cyclic Test No. 19 





Figure 3-14. Idialized Freezing/Thawing Cycle "C" 

was reestablished. This can be seen in Figure 3-15 which shows that 
depriming occurred apparently when the gas/vapor front reached the pre- 
viously frozen condenser section. 

In the last three cycles of the test program, the heat pipe, under a 
continuous load of 130 watts, was subjected to condenser warmup rates 
ranging from 0.2C/min to 0.4C/min. Cycle 25 with the higher CWR of 0.4C/ 
min was normal. Burnouts were observed in the next two cycles. The last 
heat pipe failure during cycle 27, however, is worth noting since it 
occurred 165 minutes into the warmup period when the CWR was increased 
from 0.2C/min to 5C/min. The temperature history of this cycle is shown 
in Figure 3-16. Although the cooling period of the cycle is not shown in 
this figure. Figure 3-2 indicates the ice front prior to warmup was located 
somewhere between sensor 9 and 10. Figure 3-16 shows the gas/vapor front 
at the end of the slow warmup period was located between sensors 8 and 9. 

As the warmup rate was increased, it can be seen that the gas front 
rapidly moved between sensors 9 and 10 seemingly reaching the previously 
frozen section and triggering a burnout. 
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Figure 3-15. Temperature Histories During SN009 Heat Pipe Cyclic Test No. 23 





Following this final burnout, the heat pipe was leveled and left 
undisturbed over the weekend. Fifty-five hours later, when 130 watts were 
applied to the evaporator, the heat pipe failed to hold the load, an event 
seemingly pointing to the survival of bubbles from the last freeze/thaw 
cycle as the probable cause of artery depriming. The temperature histories 
of cycles not shown in this section can be found in Appendix B.2. 

3.4 SIGNIFICANCE OF TESTS RESULTS 

The results of the cyclic tests on the SN009 heat pipe clearly estab- 
lished that thawing the frozen condenser of an arterial heat pipe under 
high load is an operating mode that results in arterial failure with sig- 
nificant frequency. Out of twenty one freezing and thawing cycles under 
high load, sixteen cycles resulted in heat pipe failures. Although thaw- 
ing the condenser under high load enhances the probability of artery 
depriming, it is not a necessary condition. This proposition is supported 
by the fact that out of three cycles in which the condenser was thawed 
under no heat load one cycle (number 23) resulted in artery failure sub- 
sequent to thawing when the heat load was reestablished on the pipe. 

In twenty four freezing and thawing cycles no heat pipe failures 
occurred during condenser freezing. The lack of failures during freezing 
when the suction freezeout mechanism has its greatest potential to deprime 
the arteries seems to indicate that none of the observed failures during 
the test cycles can be solely attributed to this mechanism. Even though 
bubbles are known to be released during freezing they are immobilized 
inside the ice structure rendering a failure during freezing due to bubbles 
an unlikely event. 

All the observed anomalies occurred during or following condenser 
warmup periods. The fact that the inactive condenser was either completely 
thawed or still frozen along a considerable length when deprimings occurred, 
it can be argued that the freezing blowby mechanism could not account for 
any of the heat pipe failures observed during the cyclic tests. 

Since neither suction freezeout nor freezing blowby were the sole 
cause of any of the observed heat pipe failures, it can be concluded only 
that all the anomalies during the cyclic tests on the SN009 heat pipe were 
the result of the freezing/thawing bubble mechanism. 
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The freezing and thawing process is required for bubbles to form 
inside the arteries. Periodic freezing and thawing replenishes the bub- 
ble population which enhances the probability of artery depriming due to 
bubbles. The results of three normal successive cycles without freezing 
supports the above. 

The presence of bubbles in the arteries is necessary for depriming 
to occur. However, depriming does not necessarily occur when bubbles are 
present. This statement is supported by the fact that out of 24 freeze/ 
thaw test cycles seven were normal (i.e., no failure occurred). These 
test results indicate that the occurrence of the observed heat pipe fail- 
ures is random, which is consistent with the demonstrated statistical 
nature of the bubble mechanism. 

During the freeze/thaw process bubbles are generated in the gas- 
blocked condenser section. It was postulated at the outset of the cyclic 
tests that the active condenser length had to expand into the previously 
frozen section. However, in roughly half the failures the active conden- 
ser length had not yet reached the previously frozen section. 

The results of the tests seem to indicate that there exists some cor- 
relation between the condenser warm-up rate and the incidence of heat pipe 
failures. The reduction of the condenser warm-up rate used in the first 
five cycles by approximately 50 percent was followed by the only two con- 
secutive normal cycles observed during the tests. In addition, during 
cycle 27 depriming occurred almost two hours into the warm-up period fol- 
lowing a rapid and large increase in the condenser warm-up rate. 

Deprimings observed during cycles with high warm up rates occurred 
after the condenser had completely thawed and soon after the active con- 
denser length reached the previously frozen section. On the other hand, 
most heat pipe failures during cycles with low warm up rates occurred 
while portions of the condenser were still frozen and before the active 
condenser length ever reached the previously frozen section. The latter 
failures are attributed to movement of bubbles into the active condenser 
helped by suction freezeout and/or Marangoni flow. 

The results of several unsuccessful attempts to reprime the heat pipe 
after the first freeze/thaw cycle indicate that a depriming does not clear 
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out all the bubbles. Consequently bubbles form during a single freeze/ 
thaw process can result in repeated successive failures due to residual 
bubbles. 

Numerous deprimings observed during the tests occurred 10 to 30 min- 
utes after freezing conditions in the condenser ceased to exist. Follow- 
ing the last freezing and thawing test cycle, the SN009 heat pipe deprimed 
52 hours later when the normal heat load was reestablished. This heat pipe 
failure can only be attributed to be the result of bubbles generated during 
the last freeze/thaw cycle of the test program. 

These results indicate that ice-generated bubbles can last a long 
time, and therefore freezing at the very time of the CTS anomalies was not 
necessary. Thus the anomaly on day 253 can be attributed to ice-generated 
bubbles from day 252 or earlier. 
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4.0 CONCEPTS TO AVERT ARTERY DEPRIM IN6 


4.1 NONCONDENSIBLE CONTROL GAS SELECTION 

Unquestionably the simplest and most effective concept to avert arte- 
rial depriming is the selection of a control gas that would preclude the 
formation of bubbles when the liquid methanol undergoes freezing and thaw- 
ing. This concept is unfortunately unrealizable, for the freeze/thaw pro- 
cess, being a degasifying process, will result always in the formation of 
bubbles regardless of the noncondensible control gas selected. 

Nucleation experiments performed during Phase I, showed that a large 
number of minute bubbles are initially released during freezing and thaw- 
ing. Most of them rapidly disappear owing to their size, but others coalesce 
forming fewer but larger bubbles whose lifetimes strongly depend on their 
sizes. Larger bubbles have longer lifetimes which enhance the probability 
of further coalescence which results in bubbles of even longer lifetimes 
which can then survive until the gas/vapor front moves into previously 
frozen sections or until suction freezeout and/or Marangoni flow succeed 
in making them grow and transporting them into the active section of the 
heat pipe. 

The above observations indicate that reducing bubble lifetimes by 
proper selection of a control gas can decrease to a certain extent the 
probability of bubble coalescence and of their induction into the active 
section of the heat pipe. The results of the bubble lifetime analysis in 
Section 2.0, show, for example, that a spherical CTS control gas bubble 
(90% N2/10% He) of 0.8 mm in diameter (i.e., half the CTS artery diameter) 
under some assumed conditions has a lifetime of about 19 hours compared to 
a lifetime of less than one hour for a methane bubble of the same size. 

The rapid dissolution of methane bubbles suggests the selection of methane 
for the control gas. 

Although the selection of a new control gas, such as methane, will 
result in bubble lifetimes significantly shorter than in the CTS heat 
pipes which will diminish to a certain extent the probability of deprim- 
ing, this concept cannot guarantee that depriming will be averted. However, 
the use of methane in conjunction with other steps to avoid depriming 
should enhance the effectiveness of those steps. 


4.2 OPERATIONAL PROCEDURE FOR START-UP OF A FROZEN CONDENSER 


Since the bubbles leading to arterial depriming are formed by a 
freeze/thaw process, the avoidance of freezing will avert artery deprim- 
ing. Selection of sun angle and radiator radiation properties or the use 
of heaters could prevent freezing in many missions. However, in missions 
to outer planets, avoidance of freezing may be costly or even impossible. 
Even when the condenser has been frozen, it is possible to thaw the heat 
pipe and restore it to a primed and serviceable state, provided proper 
procedure is followed. 

The key to successful restoration of a previously frozen heat pipe is 
dissolution or venting of the freeze/thaw generated bubbles. Venting is 
discussed in Section 4.3„ Dissolution requires leaving the condenser in a 
quiescent state somewhat above the freezing point for a time long enough 
to allow the bubbles to collapse. This time depends upon the control gas, 
the size of the largest bubble, and the pressure and temperature in the 
gas-blocked condenser. The advantage of methane as a control gas has been 
discussed. It has a high solubility-diffusivity product leading to more 
rapid bubble dissolution.. Since the bubble agglomeration mechanism appears 
to be statistical, and it is not clear that a laboratory test at 1-g is a 
good indicator of agglomeration at zero-g, it is presently difficult to 
predict with certainty the size of the largest bubble and thus the time 
required for it to collapse. However, start-up can be attempted earlier 
than necessary merely to test whether a sufficiently long time has been 

r 

allowed to elapse. If depriming occurs, the heat pipe is returned to a 
quiescent state for more time to pass. In the attempt at operation, at 
least one bubble was vented. Thus successive attempts at start up contrib- 
ute to the bubble clearing process. 

High pressure during the bubble collapse period hastens it. With an 
actively controlled heat pipe, the gas reservoir and evaporator could both 
be heated somewhat to raise the gas-blocked condenser pressure. It should 
be noted that the restart procedure would have to be used each time a seg- 
ment of the gas-blocked condenser freezes and is thawed. 


4.3 MECHANICAL MODIFICATIONS 


In a gas-controlled arterial heat pipe it has been found that a prim- 
ing foil at the end of the evaporator is necessary for successful priming 
[5]. Without it, gas in the pipe can be trapped inside the artery during 
priming. With the priming foil and with a gentle heating of the evaporator 
bubbles trapped in the evaporator are transported to the foil and vented. 

The priming foil shown in Figure 4-1 consists of a very thin metallic 
membrane with one or more holes no larger than the effective pumping pore 
diameter needed for the capillary flow in the active pipe. The size of 
the vent hole and the thickness of the foil must be such that the meniscus 
of a bubble of arterial size contacts the meniscus attached to the outer 
corners of the vent hole. Meniscus coalescence then ruptures the liquid 


i 
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Fiugre 4-1. Segment of an Artery with a 
Priming or Venting Foil 
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film bridging the vent hole and allows the gas to be vented. It should be 
noted that the heat load during venting must be less than the open-artery 
capacity of the heat pipe. 

The same principle can be applied to the condenser region or even to 
the entire pipe [5] as long as the local stress does not exceed the capil- 
lary pumping capability of the artery at that point and time. Thus a loca- 
tion in the condenser (or perhaps even in the adiabatic section) can be 
identified where under restart heat load the local stress is below this 
limit. At that location a screen across the artery can be placed to fil- 
ter out and collect small bubbles tending to move toward the evaporator. 

On the condenser side of the screen a priming foil would be placed (venting 
holes located one artery radius from the screen) to allow venting of the 
collected gas should the gas volume reach arterial size. 

The bubble-filter/venting-foil concept should allow a normal starx up 
from a frozen condenser state and permit continuous normal operation of 
the heat pipe even though the gas-blocked condenser is undergoing 
periodic freezing and thawing. 
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5.0 CONCLUSIONS 

Evidence strongly suggests that bubbles generated through freezing 
and thawing of methanol in the gas-blocked condenser arteries were respon- 
sible for the CIS thermal anomalies. Out of the 24 freeze/ thaw cyclic 
tests with the SN009 heat pipe, 17 deprimings were observed. While some 
of the tests were conducted under conditions wnere probability of deprim- 
ing was high, others were conducted where the prt^bability was low, thus 
emphasizing a statistical variation in the freeze-thaw bubble mechanism. 

Most often depriming occurred while the gas-blocked condenser was 
still frozen along a portion of its length but some occurred 10 to 30 min- 
utes after thawing was completed. In one test depriming occurred 52 hours 
after thawing. 

Analytical -numerical studies of bubble collapse show bubble lifetimes 
to be strongly dependent upon bubble size and control gas composition. 

The helium in the CTS control gas mixture causes long bubble lifetimes 
when the bubbles are of arterial size. Selection of methane as a control 
gas would significantly shorten bubble lifetime for equal sized bubbles. 

At the present level of understanding, selection of methane for control 
gas cannot ensure avoidance of bubble depriming. Prevention of freezing in 
the gas-blocked condenser by active or passive control should avert deprim- 
ing. Should avoidance of freezing be inconvenient or not feasible, either 
installation of a bubble-filter/priming foil element at the condenser- 
adiabatic section interface or adoption of an appropriate thawing and bubble 
clearing procedure would probably allow successful heat pipe operation. 
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APPENDIX A.l 

EXCESS SKEW IN LOAD PARTITIONING 


The three heat pipes of the CTS VCHP system are designed to turn on 
sequentially so as to balance the load between them at full-on design con- 
ditions. However, a build up of tolerances on gas inventory during manu- 
facture, and very low sink temperature during equinox conditions may result 
in the load on heat pipe No. 1 exceeding its capacity before Nos. 2 and 3 
take up a significant share. Failure of No. 1 by the above mechanism or 
by any other could then result in a rapid transfer of load to No. 2, etc., 
perhaps yielding a "domino effect" failure of the whole system. 

The first step in investigating this was to further review the flight 
data to search for regularity in the depriming sequence associated with the 
anomalies. Toward this end, NASA LeRC personnel provided computer plots of 
pertinent flight data for subsequent review. 

Experimentally, attempts were made to deprime the SN009 heat pipe 
with rapid load increases, simulating rapid sequential load transfer as 
pipes deprime. 

The tasks performed in this investigation are summarized below. 

Flight Data Analysis to Determine Depriming Sequence 

Computer plots and tabulations of CTS telemetered data were provided 
for days 75, 82, 89, 101 and 253 of 1977. Day 39 was a so-called "normal" 
day whereas the others were the four anomaly days. 

Flight data had already been examined in a previous study effort. 
However, the new data included overlay plots of all pertinent data on 
single pages, and value'*> of aT 2_2 and aT 2 _^ in addition to provided 

earlier. It was hoped that this data format might shed additional light 
on the system behavior, particularly the sequence in which the heat pipes 
deprimed on the anomaly days. 

Unfortunately, the telemetry increment in the temperatures wasl.SC 
to 1.9C, which is of the same magnitude as the normal values of aT 
between pipes as well as the changes in temperature which accompany 
depriming. Consequently, it is very difficult to distinguish between 


i 

actual physical phenomena and telemetry noise. In view of this limitation, 
the data tentatively suggests the following; 

Day 75: It appears that heat pipe No„ 1 deprimed first, yielding 

increased values for ^^3-1* ^^^eral hours later 

heat pipe No. 2 deprimed, followed almost immediately by 
heat pipe No. 3 and the sudden rise in body temperature. 

Day 82: The data is particularly unclear this day, but as with day 

75, it appears that heat pipe No. 1 failed first, followed 
by Nos. 2 and 3. 

Day 89: No anomaly occurred on this day. However, the data sug- 

gests that heat pipe No. 2 was deprimed, causing No. 3 to 
carry a greater load. 

Day 101: The data indicates the failure sequence was 1:2:3, as with 
day 75. 

Day 253: In this case, heat pipe No. 3 did not even “turn-on" until 

the anomaly began, and apparently could not prevent it. 

Thus, it had to be deprimed from the start. In fact, the 
data suggests that heat pipes 2 and 3 deprimed first, with 
the anomalous increase in body temperature corresponding to 
the depriming of No. 1. 

Specific interpretation of the depriming sequences is covered in the 
SN009 heat pipe test discussion that follows. However, one item of note 
bears discussion. If bubble nucleation within the arteries is a cause of 
the anomalies, then the statistical nature of the mechanism should yield 
one heat pipe failed occurrences far more often than two or three heat 
pipe failed occurrences. Day 89, in which No. 2 appears deprimed, while 
Nos. 1 and 3 are primed, may be such a day. 

Rapid Increases in Load Tests on SN009 Heat Pipe 

Extensive testing was performed witn the SN009 heat pipe. The first 
series of tests were to determine whether there exists significant iner- 
tial effects on heat pipe capacity which could contribute to a "domino 
effect" failure of all three CTS pipes. That is, if heat pipe No. 1 



deprimed, most of the load it had carried would be suddenly transferred 
to No. 2 , etc. If the capacity of a heat pipe were substantially lower 
for the sudden application of load than for gradual application (due to 
fluid inertia), the probability of a sequential failure of all three pipes 
would be enhanced. 

In order to maximize the rate of load transfer into the heat pipe 
evaporator, attached heater blocks normally used to simulate the thermal 
mass of a heat source were removed and tape heaters were placed on the 
heat pipe saddle over a 30-cm long section of the pipe which was then 
properly insulated. Testing was done for both high (-18C) and low (-96C) 
sink temperatures with the heat pipe evaporator elevated 2.5 cm. Heat 
loads ranging from 100 to 150 watts were rapidly applied resulting in no 
heat pipe failures. The results indicated essentially no rate effect 
attributable to fluid inertia. This argued against the postulated 
"domino effect" failure of the whole VCHP system. 
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INTRODUCTION 

When a gas-controlled heat pipe's reservoir experiences an eclipse, 
the gas partial density within the reservoir rises, and the gas front re- 
treats until the evaporator temperature drops sufficiently to establish a 
new equilibrium. During the transient, the evaporation from the warm end 
may exceed the capacity of the axial wicking, and arterial depriming may 
occur. It Is desired to estimate the Instantaneous heat flow during such 
a transient for comparison to the burnout heat flow capacity. 


RESERVOIR COOLING RATE 

The reservoir temperature versus time may be known from measurement. 

It is set by a heat balance between the surrounds and the thermal capacity. 
Let the mass of the reservoir be M and Its specific heat c. Conductive 
coupling to the radiator Is KA/L, and Insulation is modelled with massless 
radiation shielding. 

-me dT - AoF [oT^ - a q"] - KA [T -T^] (1) 

dt 7 T 


Before eclipse q” is high, and the quasl-equilitrium temperature is given 
by setting the left-hand side to zero. After q~ Is reduce^T, Cools at the 
rate given by Eq. (1). 


FLOW VELOCITY INTO RESERVOIR 


When a wicked reservoir of volume V cools at a rate of -dT/dt* the 
vapor partial pressure Inside the reservoir falls. Since the total 
pressure Is fixed at the evaporator end (neglecting vapor flow pressure drop), 
the partial pressure of the gas rises. Furthermore, the gas temperature falls 
(assuming thermal equilibrium). These two effects Increase the gas partial 
density. The rate at which the partial density times the volume rises must 
come from the flow of gas from the heat pipe radiator Into the reservoir. 



tPj V) 



( 2 ) 


P - P-P^(T) 


(3) 


Differentiating Eq.(3) gives 



P-Pv 1 /-dT\ + 1 dPv 

ITT T I Ht j O BT" 

9 9 


/-dT\ * 1 dP 
V 3t/ RgT dt 


Substituting into Eq.(2) then gives 



+ 


dPvT /- dT\ + dP 1(4) 
?rj \ W dt j 


Let dPv/dT be approximated from the fclausius - Oapeyron equation 

V 0 
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Then Eq.(4) becomes 


V 


X 

A. 



+ Pv (Ah -1 

r ^ 



( 5 ) 


The total pressure P is the vapor pressure at the evaporator temperature T^. 


dP = dPv / dTv \ 
dt dTv \dt / 


“ iEl /-dTv\ 
dTv dt/ 


( 6 ) 


The evaporator temperature in turn is set by a heat balance upon it. The 
worst case occurs when the evaporator is well-coupled thermally to a large 
thermal mass. In this worst case dP/dt may fall only very slowly compared 
to the other terms. Measured data may be used to evaluate dTv/dt. A crude 
model is a one-capacity model. 

= «elec. Alec. Ap -A# (=t/ -ot/) (7) 


The last term in the equation models heat leak to surrounds at temperature 


HEAT FLOW 

The instantaneous heat flow through the pipe Q^p is fixed by the 
length of active condenser and the rate of advancement of the gas front. 

• • • . 

^hp “ ^static ^dynamic (8) 


The static load is that caused by condensation on the gas-free condenser 


wall 


^static * ” ^*^^v “ S. Sc ) 


(9) 


where the length L. Is fixed by a gas Inventory m. 

c g 

•"g ' Pg V + PggAc (L -t-c) 00) 

The dynamic load is proportional to the velocity v at which the front 
advances 

. Vamic * " (T, -V 01) 

where is the condenser fin density, c^ its specific heat, its area 
(volume per unit length), and n is not the radiating fin effectiveness, be 
cause of the fourth power dependence of the black body radiosity. The 
appropriate value of n is between the radiating fin effectiveness and the 
fourth root of that value, depending upon the ratio of 


EXAMPLE 

For example consider the CTS heat pipe 1 on Day 75. The following 
data are available 


Time 

Vapor Tei 

t 


min 


460 

27.5 

470 

28 

480 

29 


Gas -Blocked Temp. 

^c 

°C 

-93 

-94 

-95.5 


Gas Reservoir Temp. 
T 

-68 

-71 

-74 
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From these data one can estimate -dT/dt * 18 K/hr and dT^/dt 
The following vapor pressure data pertain 


4.5K/hr. 


Temperature Pressure Latent Heat 


19500 


1.15 X 10' 


T 205 M3 1.2 X 10” 

T- 180 M.5 1.2 X 10® 

For a 90% N 2 10% He mixture the mean molecular weight Is 25.6 and Rg « 
8314.3/25.6 * 325 J/kg K. The density of the gas In the gas-blocked condenser 


19500 = 0.333 kg/m' 


Other parameters needed are 


Ah 

RT T = 205 


1.2 X 10 ” = 22.51 


Ah I 1.15 X 10” 

RT T - 300 = (260) (300) = 14.74 


fhe pertinent parameters of the pipe itself are understood to be as 


follows: 


' it *3 

Reservoir Volume V 8.70 inch'* 1 .426 x 10”^ 

Pipe Vapor Area A^. 0.1245 1n^ 8.03 x 10”® m^ 
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Equation (5) then gives 
V « 1 .426 X 10 


8.03 X 10 




V = (1.78) (4.50 X 10"®) ^ [95.1 + 1.36] (18) + (958) (4 

8 X 10"® ( 1736 + 4311 > « 0.484 m/hr = 1.59 feet/hr 


I 


+ 19500 (14.74) (4.5) 

300 


I 


! 

.5) I 


Note that* in this example* the fall of the vapor pressure in the reservoir 
was a negligible consideration but the rise of vapor pressure in the 
evaporator was quite significant. 

For purposes of an estimate nAp is taken to be 0.04 inches by 12 inches, 
Thus Equation (11) gives 

^dynamic “ (.04/12) (12/12) (1.59) (300 -180) (1.8) 


Vamlc * ' 12.25 Watts 

hr 


DISCUSSION 

While a rapid decrease in reservoir temperature has the potential of 
adding a significant dynamic heat flow to the static load* the example worked 
here indicates that the decrease must be rapid, much greater than the 18K/hr 
value used in the example. 
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Intyoduction 

When a gas -controlled heat pipe's radiator experiences an eclipse, 
the heat flow in the pipe increases until the evaporator cools to a new 
equilibrium operating condition. During the transient, the evaporation 
from the warm end may exceed the capacity of the axial wicking, and 
arterial depriming may occur. It is desired to estimate the instantaneous 
heat flow during this transient for coiiqparison to the burnout heat flow 
capacity. This memo is a conpanion to one that explored the instantaneous 
heat flow after heat pipe reservoir eclipse. 

A Sjjiplified Thermal Model . 

To explore tlie magnitude of the expected effect and to show the effect 
of the major parameters, a sinplified model is proposed. In this model the 
heat loss out of the radiator is taken to be 

vihere is the so-called script - F radiant transfer factor (equal to 
emissivity for an isolated radiator) W is the width or perimeter of the 
panel per pipe (including both sides if both sides radiate) , L is the active 
length of the vapor in the radiator, T^ is the vapor temperature, Tg is the 
environmental temperature 

T, ■ C e) 


and e£ is the fin efficiency. The vapor filled length is assumed to be 
given by 


T -T 

, V o 
^tot “2T" 


C5) 


idiere is the total length, is the tum-on point, and AT is the control 
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band width. Equation (1) may be factored into the £otb 


5 - * yffv - V 

and Eq. (3) may be introduced to write 

5 ^ " '^'e^ 

vhere 




This latter parameter is insensitive to mild excursions in and 

since the absolute tenperatures appear. In contrast Eq. (4) varies rapidly 
with Ty, because (Ty-T^) and (Ty-T^) change significantly when Ty changes, 
particularly the former quantity, and (T„-T 1 changes when changes. 

Note that this simplified model neglects the thermal capacity of the 
radiator on the grounds that it is small compared to that of the heat source. 
This neglect leads to an overprediction of thermal shock on the heat pipe from 
sudden changes in environmental tenperature T^. 

The heat source is modelled as a single lunped capacity of mass m and 
specific heat c. 

dT, , . 

vdiere Tj^ is the heat source tenperature, is the net electrical power 
dissipation, and R is the theimal resistance into the heat pipe evaporator. 

It is the heat flow into the evaporator that concerns us 

^ ■ jr ffh • V 

e 

The heat source thermal capacity me is assumed to be sufficiraitly large 
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that Tj^ remains constant during the time that the radiator and gas front 
readjust. In this respect the heat transfer through the heat pipe is 
assumed quasi-steady. 


■Analysis Based Upon the Simplified Model 

The model in its bare essenti?ils consists only o.C Eqs.(4) and (7). 
The heat pipe behavior is detennined by the evaporator resistance R_, 

w 

the radiator resistance the tum-on point T^, and the control band 
width AT. The parameters R^ and R^ may be inferred from observed temper- 
atures at, say, the full-on operating point. Equation (4) shows that 



"^v.full - "^e 
^full 


+ AT 
o e 

^full 


( 8 ) 


and from Eq. (7) 


"^h.full - '^v,full 
^full 


(9) 


Hence 


\ ' ’^v.full 

^ ^ 

^ V,full ^e 


CIO) 


Equating Eqs. (4) and C7) and rearranging gives 


ATCTj^ - Ty) -r(Ty - T^DCTy - T^) 

This quadratic equation may be solved for T.^, and the result put into 
Eq.(7). Equations (8) and (9) are also introduced ior convenience 




ATq + AT/r-.^CTo - Tg)^ + 2 (AT/r) AT^^ -t-CAT/r)^ (11) 

2r(Tjj - Tg + AT) 
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where, to shorten notation, 

' 2(T^ - ip - (T„ - T^) C 123 

Note that in the form of Eq.(ll) neither nor appear explicitly, but 
only their ratio r. 


Sample Parametric Calculations 

To show the importance of the parameter r we show three sample calcula- 
tions, one with r=l, one with r=Js, and another r»l/8. We choose the following 
nominal values: 


Tq = 20'^C 

AT = S°C 


T 

e,max 

T 

e,mm 


-20°C 

-110°C 



r = 1, 1/8 


In order to keep the same nominal value of R^ as T^ is changed, one multiplies 
by the ratio of (Tq * ^T)/CTq ‘ The following results are 

obtained: 

TABLE 1 

EFFECT OF EV^APOR\TOR-TO-RADIATOR RESISTANCE RATIO 
ON HEAT PIPE OimOAD UPON RADIATOR ECLIPSE 


r 

°C 

CQ/Qpun][(To - Tg + aT)/(T 

0-1 

1 

-20 

[0.1981][45/45] = 0.1981 

, B% 


-no 

[0.07133E135/45] = 0.2140 

0.25 

-20 

[0.6074][45/45] = 0.6074 1 

, 28X 


-no 

[0.2571][135/45] = 0.7714 J 

0J25 

-20 

[0.9384][45/45] = 0.9384 

1 45« 


-no 

[0.4550][135/45] = 1.3651 ^ 
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Nonlinear Effect 

In the preceeding the effect of radiator eclipse was explored in ’ 
the linearized limit. When the sink tenperature changes greatly, the 
heat flow is affected by the nonlinear variation of in the radiation 
terms. A numerical calculation is easily made despite the nonlinearity. 

We had 

^ * ''sCR "f '"-tot (V • 


and 


Q • ^ rth - V 

e 


(14) 


with 




(15) 


Equating (13) and (14) and introducing (15) gives 


T - T 
V o 


(Th - V ■ -IT 


T - T 

V e 


^ * T^.l) 


(16) 


Equation (15) is understood to apply to both the original condition before 
eclipse when * T^. and Tg = and to the conditions after eclipse 
vdien Ty and are colder. 

Given the value of Tj^, the original sink condition the turn- 

on point Tg, the band width AT, and the value of r based ipon the 
original conditions, Eq. (16) is used to find ^ and from Eq. (14). 
Then^vrfien goes to equation is lased again to find Ty. 2 » 

and Q 2 is found from Eq.- (14). * 
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Consider the example used before: 


T - 20°C 1 “ -20°c AT - 5°C T, = 30°C 

o n 

r = 1/4, - -110°C 

• c 

In the lineal ized limiit we obtained Q 2 /Q 2 “ 1*28. Allowing for the non- 
linearity gives the following values: 


Condition 1, By Binary Search 


°C 

30 

20 

25 

22,5 

23.75 

23.125 

23.4375 

23.2813 

23.2031 

23.1641 

23.1836 

23.1738 

23,1690 

23.1665 


LHS 


0 

10 

5 

7.5 

6.25 

6.875 

6.5625 

6.7188 

6.7969 

6.8359 

6.8164 

6.8262 

6.8310 

6.8335 


FHS 

°C 

25 

0 

11.25 

5.31 

8.20 

6.738 

7.4658 

7.1008 

6.9193 

6.8287 

6.8739 

6.8513 

6.8401 

6.8344 


Ratio of Q 2 /Q 1 


Condition 2, By Binary Search 


T 

V 

ms 

RHS 

°C 

30 

0 

46.35 

20 

10 

0.00 

25 

5 

21.55 

22.5 

7.5 

10.38 

21.25 

8.75 

5.095 

21.875 

8.125 

7.714 

22.1875 

7.8125 

9.0422 

22.0313 

7.9688 

8.3767 

21.9532 

8.0469 

8.0455 

8.046/6.834 « 

1.18 
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This ratio differs appreciably frcsn unity, but somewhat less so than the 
value obtained with the linearized equation. 


Conclusion . 

The importance of the ratio of the thermal resistance between the 
equipment and the evaporator to the thermal resistance between the radiator 
and the evironment is illustrated. When the ratio is large (near one), 
the transient load is only a few percent of the base load. When the ratio 
is small, the transient load Is a significant fraction of the base load. 


A-18 


APPENDIX A. 2. 3 

CONDENSER AND/OR RESERVOIR SHADOWING TESTS ON SN009 HEAT PIPE 


SN009 heat pipe tests were directed toward examining two potential 
artery depriming modes: 1) rapid chilldown of the "ondenser, and 2) rapid 

chilldown of the reservoir. A large evaporator mass was a requisite con- 
tributor during these hypothesized depriming modes and, accordingly, four 
kilograms of aluminum were attached to the SN009 evaporator. 

Two series of rapid condenser chilldown tests were run. In the first, 
the heat in the condenser sink was cooled to -73C (-100F), while the con- 
denser was heated to maintain -18C (OF). The sink and condenser were 
coupled through 0.30 centimeters of cork over a 162 in^ (0.1045 m^) area. 
With the heat pipe operating at 100 watts (well above open artery capacity), 

the condenser heater was turned off allowing the condenser temperature to 

fall approximately 2.3C/min (5F/min). When no depriming occurred, the 
tests were repeated with 100 and 150 watts. Depriming at 160 watts under 
steady state was verified. Thus, cooling the condenser at approximately 
2.8C/min (5F/min) produced no de priming. 

In the second series, the sink was -129C (-200F). When the con- 
denser heater was turned off, a cooling rate of approximately 3 . gc/min 

(7F/min) was achieved. When no depriming occurred, the test was 

repeated with 125 watts, and again no depriming was observed. However, 
when an attempt was made to raise the power level to 150 watts, depriming 
occurred (at the steady condenser temperature of-18C (OF). The pipe 

was reprimed, held 100 watts, but failed to retain prime at 125 watts. A 
third attempt resulted in failure to reprime. It is thought that ice- 
generated bubbles may have been formed in the earlier transient cooling 
tests, and these bubbles caused the repeated deprimings. Depriming due 
to rapid chilldown of the condenser was not demonstrated. 

Rapid chilldown of the reservoir was achieved by blowing vapors from 
boiling liquid nitrogen through a cooling coil block attached to the 
reservoir. With 100 watts power at the evaporator and the condenser at 
-13C (OF), the reservoir was cooled repeatedly at successively higher 
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rates. Depriming did occur at a rate of -3.2C/min (-5.6F/min). Deprim- 
ing did not occur at rates of -0.8, -1.2, -1.6 and -2.5C/min (-1.5, -2.1, 
-2.8 and -4.5F/min, respectively). 

Although depriming due to reservoir chilldown was demonstrated, the 
rates required were substantially higher than those indicated by flight 
data. This argues against reservoir and/or condenser shadowing as the 
cause of the CTS anomalies. 


APPENDIX A. 3 

FREEZING BLOWBY TESTS ON SN009 HEAT PIPE 


Freezing blow~by was hypothesized as a possible depriming mechanism. 
This requires the formation of an incompletely frozen slug bridging the 
vapor core during a transient in which the pressure on the evaporator side 
of the slug is increasing with respect to the reservoir side. The result- 
ing pressure difference across the slug can blow liquid from the evaporator 
to the reservoir side and deplete the evaporator inventory. 

To explore this mechanism, it was first necessary to determine how to 
perform a meaningful 1-g test in the absence of the natural slugging of 
excess liquid which occurs in 0-g. A successful technique to form an ice 
plug in heat pipe SN009 was developed. The procedure is as follows: 

a) Instrument the heat pipe with temperature sensors along its 
length to able monitoring the location of the gas front. 

b) Apply sufficient power to turn the heat pipe on i.e., to 
move gas front into the condenser section. Maintain the 
sink temperature a few degrees above the methanol freez- 
ing point. 

c) Midway along the inactive portion of the condenser, locally 
code a short ('V/2.5 cm) section of the pipe to subfreezing 
temperatures. Periodically raise the reservoir-end of 
pipe to force excess liquid to pass over the frozen pipe 
section in order to enhance ice build-up. 

d) Apply periodically, a power pulse to gas reservoir to induce 
a transient temperature/pressure increase in the reservoir 
while monitoring the temperature profile about the location 
of the gas front. Formations of all ice plug will decouple 
the reservoir from the evaporator side of the heat pipe. 

The presence of an ice barrier blocking the arteries and 
vapor spaces can then be determined from the temperature 
profile along the evaporator side of the plug which should 
not respond to an induced temperature change in the gas 
reservoir. 

Experiments were then directed toward investigating the hypothesized 
freezing blow-by depriming mechanism. The SN009 heat pipe test was modified 
such that a small, independently chilled, ice plug heat sink was attached 
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to the condenser midway along its length. The basic test approach was as 
follows: 

a) The ice plug heat sink was employed to generate a solid plug 
midway along the condenser. 

b) The evaporator power was raised to a value in excess of the 
open artery heat pipe capacity, raising the temperature and 
pressure on the upstream side of the ice plug. 

c) The ice plug was thawed, allowing the pressure differential 
to relieve itself by blowing liquid and gas through the 
thawed region. 

If there is liquid communication between the region to thaw first and 
the arteries, the above sequence of events was hypothesized to lead to 
depriming of the arteries by pumping them dry. The experiment was run so 
that this liquid communication condition existed and, in fact, the arteries 
did deprime. The experiment was run three times, always with the same 
result. Thus, it has been clearly established that freezing blow-by is a 
bonafide depriming mechanism and a candidate to explain at least some of 
the observed anomalies. 

To enable modelling the blow-by mechanism, the variable conductance 
heat pipe subroutine was expanded to include the effects of an ice plug 
forming in the condenser section of the heat pipe and open or deprimed 
artery performance. Such a plug decouples the reservoir from the active 
portion of the heat pipe and alters its control characteristics. It also 
provides the basis for the freezing blow-by depriming hypothesis; i.e., a 
pressure difference is established across the ice plug which blows the 
liquid through as the plug begins to thaw resulting in artery depriming 
which causes the heat pipe to continue operating with reduced or open 
artery capacity. 

In Appendix A. 4. 2 a transient thermal analysis of a simple VCHP/ 
radiator system model is described which shows the effects of an ice plug 
formed in the heat pipe condenser and subsequent blow-by depriming of the 
arteries on the performance of the sample VCHP system. 


APPENDIX A. 4.1 

BRIEF REVIEW OF LeRC THERMAL STUDIES 


In support of Contract NAS3-21740 (Study of Liquid Dynamics in CTS- 
Type Heat Pipes), NASA LeRC performed thermal analyses on a 'lump' param- 
eter thermal model of the Transmitter Experiment Package. These efforts 
we^’e conducted by Louis Gedeon of LeRC and are described in two NASA docu- 
ments entitled "Comparison of Predicted Versus Measured Temperatures for 
CTS Thermal Anomalies," dated October 1, 1979, and "Modifications to CTS 
Thermal Model for TRW Contract NAS3-21740," dated January 9, 1980. A 
review of this analytical task is presented in the sequel. 

The CTS model includes the VCHP system, the traveling wave tube, a 
power processing system baseplate, the spacecraft forward and part of the 
south panels, the antenna covers and skirts, and the solar array pallet. 
TRW's Systems Improved Numerical Differential Analyzer (SINDA) computer 
program and a Variable Conductance Heat Pipe model subroutine (VCHP2) were 
utilized to solve the CTS model. 

During the first transient calculations with a modified VCHP system 
model, it was found that the heat pipe model subroutine (VCHP2) was not 
suitable for transient calculations. An analysis at TRW of the solution 
algorithm in VCHP2 revealed that the solution scheme, in view of the non- 
linear nature of the model, was inadequate since it was susceptible to 
becoming numerically unstable. As a result of this analysis, a new solu- 
tion method was formulated which proved to be successful in circumventing 
numerical instabilities. 

This solution approach was incorporated into a new version jf the vari- 
able conductance heat pipe model subroutine renamed VCHPDA. In addition, 
two new features were introduced into this version which permit simulating 
the effects of an ice plug formation in the heat pipes and open or deprimed 
heat pipe capacity on the performance of the whole VCHP system. The 
description of the ice plug and open artery analytical models and the solu- 
tion scheme, as well as, a listing of VCHPDA written in FORTRAN IV computer 
language are presented in Appendix A. 4. 2. 
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Subsequent transient runs of the CIS model using »,i*e new subroutine 
VCHPDA were normal, except that the predicted temperaiuresJ were lower than 
flight temperature data. In order to improve the temperature predictions, 
additional modifications of input parameters were introduced into the CIS 
model. Among the changes made were: 

1) The solar absorptance and total hemispherical emissivity on 
the VCHP system radiator were increased and decreased 
respectively. 

2) The reflected sun heat load on the VCHPs radiator was 
assumed diffuse and was multiplied by a factor varying 
with time. 

3) The south panel emissivity was increased. 

4) Numerous view factors and shadowing parameters were 
modified. 

Although some of the above ch'-figes are deemed somewhat physically 
unrealistic, they resulted in general, in temperature predictions remark- 
ably close to flight data. Transient runs were made for normal day 89 and 
for periods proceeding and including the anomaly of days 75, 82, 101, and 
253. The initial conditions were those obtained from steady state calcula- 
tions using the conditions prevailing prior to the start of a spacecraft 
eclipse. 

For each anomaly day a minimum of two transient runs were made, one 

assuming normal heat pipe operation prior and during the period of the 

anomaly, with considerations for ice plug formation, the other run assum- 
ing open artery capacity conditions prior and during the anomaly. Addi- 

tional calculations were made for day 75 and day 253. For day 75, lower 
open artery capacity values were assumed, and for day 253 a third run was 
j'serformed in which normal heat pipe operation was assumed from eclipse 
until the observed time of the anomaly at which time the heat pipes were 
set for open artery or deprimed capacity operation. 

Some significant results from these transient analyses were: 

1) On day 253 freezing was not predicted at the time of the 
anomaly. 

2) The maximum calculated heat pipe loads never exceeded 
80 watts. 
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3) The coldest calculated temperature for heat pipe number 1 
always occy 'red near the reservoir end of the condenser 
where temperature sensor HPT5 is located (see Figure 1-3 
in text). Temperatures for heat pipes HP2 and HP3 at 
similar condenser locations were predicted to be a frw 
degrees colder. 

4) Formation of ice plugs in the three heat pipes were calcu- 
lated on day 82 prior to the anomaly. 

5) The profile of the tube body temperature excursion is 
closely matched on day 253 when the heat pipes were set 
from normal to open artery capacity operation at the 
observed time of the anomaly. 
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APPENDIX A. 4.2 


GENERALIZED VARIABLE CONDUCTANCE HEAT PIPE MODELING 


This appendix presents analytical models, sample problem solution, and 
listing of subroutine VCHPDA. 
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I 

f As the result of efforts carried out under the 1973 IR&D 

L 

i Advanced Thermal Control program, a method to analytically simu- 

I late a gas loaded, variable conductance heat pipe (VCHP) was 

i developed and is documented in Reference 1. 

. A FORi’RAN IV computer code was implemented to numerically 

solve the analytical model. This code which is suitable to 
' interface with general thermal models involving heat pipes, has 

been used successfully in steady state thermal simulation of 
!' systems, e.g., the transmitter experiment package (TEP) variable 

conductance heat pipe system (VCHP) in the Canadian Technology 
Satellite (CTS). Subsequent attempts, however, to utilize the 
VCHP2 subroutine for transient thermal analyses yielded unsuc- 
cessful results due to numerical instabilities arising during 
the gas-vapor front location/vapor temperature calculations. 

Recently, as part of the efforts to investigate the thermal 
anomalies in the CTS spacecraft, the solution scheme in VCHP2 
was analyzed in order to uncover the source of the numerical 
problems. The analysis revealed that the iterative approach 
of the subroutine to obtain a converged solution starting from 
an initial guessed gas length/vapor temperature, was inadequate 
in view of the non-linear natureof the model. 

As the result of this analysis, a new analytical method 
was formulated which under most physically realizable operating 
modes, unconditionally circumvents numerical instabilities, 

i 
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In addition, two new features were introduced into the VCHP 
model which permits simulating the effects of an ice plug for- 
mation and open artery capacity on the performance of a VCHP 
system. The appropriate FORTRAN IV logic to solve the analytical 
model was incorporated into a revised version of the VCHP2 sub- 
routine, now identified as VCHPDA. A logic flow chart and listing 
of the computer program are attached. 


Basic Analysical Model; 


Referring to the configuration of a typical, wicked-reservoir, 
VCHP in Figure 1, and invoking the ideal gas law and the "flat 
front" gas theory, the cumulative distribution of the number of 
moles of non-condensible gas along the length of a VCHP can be 
written as: 


Ng(z) = (Pj(z) 




Pj(z) - Py(z) 
■R^tTY] 


A(z) 


dz 


( 1 ) 


for z - Lg 
where 

Lg - location of gas-vapor front or "gas length" 

Vp - reservoir volume 

T|^ and T(z) - temperatures at vapor-liquid 
interface in reservoir and along pipe, respectively. 

R - universal gas constant 

A(z) - gas/vapor space cross sectional area 
P D P,(z) - partial vapor pressure of working 

V 5 K V 

fluid at Tj^ and T respectively. 

P.^(z) - total pressure in the VCHP defined in Eq . (2). 


Pf(z) = Pv(Tv(z)) 


( 2 ) 


where T (z) is the vapor temperature in the active section of 
the pipe. 
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I - I 

T^(z) = f h(z') T(z') P(z') dz/f h(z') P(z') dz' (3) 
z z 

where h(z) is the local heat transfer coefficient and P(z), the 
wetted perimeter. Equations 1, 2 and 3 form a complete non-linear 
coupled set of equations which enable the calculation of the 
length of the gas-blocked region and the vapor temperature in 
the active portion given the temperature distribution along the 
wall of a VCHP whose configuration conforms to the one sketched 
in Figure 1 . 


Solution Scheme: 

In order to illustrate the solution procedure incorporated 
into the current version of subroutine VCHP DA, it is convenient 
to rearrange Equation (1) as follows: 

Ng(z) = Pj(z) F(z) - N^(z) (4) 

where 

R^T(z' ) * 


Nv { z ) 





Py(z' ) • A(z' ) 
R • T(z“^ 


• dz ' 


( 6 ) 


Resorting to a numerical method to evaluate the integrals F(z), 
Nv(z) and Pj(z) are calculated stepping out from the reservoir 
in Az increments. At each step, Ng(z) is evaluated and compared 
to the total gas inventory, N^. The above procedure is continued 
until Ng(z) ^ Nj at which point a quadratic binary search scheme 
is applied to find the root of the transcendental equation 
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(7) 


Ng(z) - = 0 

which is satisfied when z = Lg. 

To initiate this iterative process, a gas length is approxi- 
mated by linear interpolation, i.e., 

Ng(z+Az) - Ny 

'-9 ' ^ * NiT2+Szl"'-“N9T2 !' ■ 


The search is terminated when a prescribed convergence criterion 
is met, which in the current VCHPDA version is defaulted to 

I f-'g(z) - Nt I . 


Solution to Equation (7) allows the performance of the partially 
gas-blocked VCHP to be quantified. 

The thermal conductance between the vapor and the pipe 
walls, G(z), can be calculated. 


G(z) = U(z-Lg) h(z) p(z) 
where [J is a Heavyside-type function defined as 


u 


0.0 for (z-Lg) 0.0 
1.0 for (z- Lg ) - 0.0 


Q = 


where H is similarly defined as 

0.0 for I Ty(Lg) 


H = 


1.0 for 


{ T„(Lg) - 


T(z) I <0.0 
T(z) } ^ 


0.0 


( 8 ) 


and the heat load on the VCHP 

J G(z') • (T^ag) - T(Z')) • H(Ty{Lg) - T(z')) • 
1-g 


dz 


(9) 
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The heat load on the heat pipe in terms of power-length can 
readily be calculated 


where 


QL = Q • 


L 


eff 


^eff 


(L -Lg) 

— ^ + L 



( 10 ) 


( 11 ) 


The above solution scheme holds provided the VCHP is under load. 
In the event that Ng(z) < Nj for z - Ly, the gas length Lg = Lj 
and the total pressure in the VCHP is calculated by either 
Equation (1) or (4). Using the latter Equation, the pressure 
P.J. is given by 


Pf = { N' j - Ny(L^) I / F(L.^) (12) 


Ice Plug Modelling 

Owing to excess fluid inventory, usually present in a VCHP, 
operation of the heat pipe under subfreezing sink condictions 
can result in the formation of an ice barrier in the condenser 
section effectively decoupling the active portion of the pipe 
from the gas reservoir. Formation of such an "ice plug" can 
have a significant effect on the performance characteristics of 
a VCHP. Furthermore, it can, under certain operating conditions, 
lead to depriming of the arteries by what is referred to as the 
"freezing blow-by" mechanism. In order to enable to simulate 
the consequences of such an event, the ice plug model has been 
developed and is incorporated as an optional operating VCHP mode 
in subroutine. 

It is assumed that in zero-gravity the bulk of the excess 
fluid resides as a liquid mass in the gas reservoir and as a 
slug blocking the vapor spaces in the pipe section at and near 
the coupling between pipe and reservoir. During a transient 
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analysis, formation of an ice plug in the condenser is subject 
to the following constraints: a) the gas reservoir temperature 
must be above freezing conditions, b) the history of the VCHP 
must show a net decrease of the gas inventory in the reservoir, 
i.e., gas front and, presumably liquid slug, movement toward 
the evaporator, and c) subfreezing temperatures at some point 
along the condenser and/or adiabatic sections. 

These conditions for ice plug formation can be stated 
mathematical ly , 

Tr > Tp p {Freezing temperature of working fluid) (14) 

8Ngj^/9t < 0 (15) 

T(z) < Tp p_ ( 16 ) 

Provided conditions given by Equations (14) and (15) are met, 
the solution procedure subsequently involves stepping out from 
the reservoir end of the pipe in Az increments in the search of 
subfreezing temperatures. The point at which conditions. 

Equation (16), is first satisfied establishes the location of 
the ice plug which is currently modelled as a barrier of negli- 
gible thickness. The location of the ice plug, L , is defined by 

r 

Lp = z where T(z) - Tp p and T(z-dz) > Tp p (17) 

Concurrent calculation of the cummulative gas inventory up to 
z = L (Equation (1)) enables to determine the distribution of 

r 

the gas on the reservoir and evaporator side of the ice plug 
which respectively are: 


Ngp = Ng(Lp) 


(18) 


Ngp = Hj - Ngp 


(19) 
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These parameters are stored and used at subsequent times to 
analyze the characteristics of the two decoupled regions in 
the VCHP. The solution procedure is then as follows: 

I. Calculate total pressure on reservoir side of the 
plug using Equation (12) 


P 


T,R 


NQr - Nv(Lp) 


/F(Lp) 


( 20 ) 


II. Perform analysis on the evaporator side of the ice 

plug resorting to Equations (1) through (12) and 

appropriately changing the limits of integration 

and reservoir volume in order to account for the 

decoupling of the two regions, i.e., set Vd = 0 

< ^ 

and evaluate the integrals fron Lp to z - Ly. 

In addition, calculation of the effective length by Equation 

(11) requires the condenser length, L , be reduced by L . 

c p 


L 


I 

c 



( 21 ) 


The above procedure yields an estimate of a significant new 
parameter: the pressure differential across the ice plug, which 
can result in depriming of an arterial VCHP when liquid contact 
is established between the two regions. 


Open-Artery Capacity Modelling 

If a mechanism such as the one described above leads to 
depriming of an arterial VCHP, the designed capacity of the 
heat pipe will deminish to what is usually referred to as "open- 
artery" capacity. In such an operating mode, capillary pumping 
pressure is substantially reduced and continuous operation of a 
VCHP at high loads after depriming will generally result in 
partial evaporator dry-out. 

In order to simulate the performance of a VCHP with 
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open-artery capacity, a model was developed and is currently 
incorporated in subroutine VCHPDA as an operational operating 
mode. To exercise this option, a new constraint is introduced 
into the overall VCHP model. The power-length capacity of the 
heat pipe, defined in Equation (10), cannot now exceed a pre- 
scribed capacity (QL)^, i.e., 

QL(z) - (QL)^ (22) 


The effective length, is calculated by 


L 


eff 



+ L 


e 


(23) 


where, L'^, is the evaporator length reduced by the length 
of the dry section of the pipe, 


L'n = L, - L , 
e e d 


Lh = L-r - L,, 
d T w 


(24) 

(25) 


where, L^, is the length of the wet portion of the pipe measured 
from the reservoir end. 

Rearranging Equation (22) as 


QL(z) - (QL)^ = 0.0 


(26) 


the solution procedure involves finding the root of the above 
Equation which establishes the length of the wet section of the 
pipe and hence, the dry portion of the evaporator. 

The general procedure which requires a double- iterative 
scheme is as follows: 

I. Initialize = Lj (i.e., = 0.0) and Lg = Lj 

II. Calculate gas length, Lg, and vapor temperature, T , 
such that Equation (7) is satisfied 
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III. Calculate QL(z) for z = L and check whether constraint 
(Equation (22)) is met. if the cneck is positive, the 
solution to the model has been attained and the procedure 
is terminated; otherwise, perform step IV. 

IV. Holding Lg fixed, and stepping out from the evaporator 
end of the pipe in (-Az) increments, recalculate the 
vapor temperature using 

z z 

Tv(z) = f h(z') T(z') P(z') dzV f h(z-) P(z') dz' (27) 

Lg Lg 

and subsequently calculate QL(z). 


Repeat steps III and IV until QL(z) - (QL)^ at which point 
a quadratic binary search is performed to find the root of Equation 
(26). To initiate the iterative process, the wet length, L^, 
is approximated by linear interpolations as follows: 


L 


w 


= z 


(QL)^ - QL(z) 
QL(z+Az) - QL(z) 


Az 


(28) 


The search is terminated when a prescribed convergence 
criterion is satisfied which is defaulted in VCHPDA to 


QL(z) - (QL)^ 

“TQlT; 


- 10 


-4 


With the wet length, = z, held fixed a new gas length and 
vapor temperature are calculated confining the calculations to 
the wet section of the pipe. Steps II through IV are repeated 
until 


,k+l 


,K+1 


i 10-'' 


1 / 

where W is either the gas length, Lg, or wet length, L^, and 
superscript 'k' is the iteration number. 
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Having established the domain of the active section of tiie 
pipe, i.e., Lg - z - L^, the local thermal conductance, G(z), 
between the vapor and adjacent pipe walls can now be redefined 
to reflect the reduction of effective heat transfer area resulting 
from gas blockage and wall dry-out. 

6(z) = H(z) h(z) P(z) (29) 

where H(z) is defined as 

I 0.0 for (z-Lg) < 0.0 or (L,,-z) < 0.0 
H = . ” . (30) 

1 1.0 for (z-Lg) - 0 and (L^-z) - 0.0 


VCHPDA: Computer Program Subroutine Usage 

This subroutine solves numerically the analytical VCHP 
model described in preceeding sections of this report. Its 
solution logic is written in FORTRAN IV language which is com- 
patible with current Control Data Corporation (CDC) compilers. 

As a program subroutine, VCHPDA is meant to interact with 
general lumped parameter thermal systems. In its present form, 
however, VCHPDA is only suitable for usage in conjunction with 
TRW s Systems Improved Numerical Differencing Analyzer (SINDA). 

The subroutine must be called in VARIABLES 1. In the 
process of generating a thermal model which involves a VCHP sub- 
system, the following convention must be followed. 

a) Heat pipe wall nodes are numbered and input sequentially 
stepping out from reservoir end. 

b) Wall to vapor conductors are numbered and input sequen- 
tially as in a ) . 

c) The vapor of each heat pipe must be declared a boundary 
node . 
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Figure 2. VCHPDA Solution Flowchart 
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The calling sequence is: 

VCHPDA (Al, TN, A3, GN, A2 , TJ, TK, QN, K1 , K2, TL, A4) 
where : 

Al = Working fluid saturation pressure vs temperature array 

TN = First node in wall temperature array 

A3 = Array of vapor to vail conductors values with no 
gas input in order starting from the reservoir 
end (Btu/hr-°F) 

GN = First conductor in vapor to wall conductor array 

A2 = Heat pipe characteristics array (see below) 

TJ = Reservoir temperature node 

TK = Control temperature node (for a heated reservoir 
system ) 

QN = Reservoir heat input identification (for a heated 
reservoi r system ) . 

K1 = Constant set to zero, fixed point 

K2 = Usage flag; 1 = normal calculation usage 

-1 = printout of heat pipe data only, 
used in output calls 

TL = Vapor temperature node 

A4 = Array of node lengths starting from reservoir 
end (inches) 


Inputs for 
A2(l ) 
;\ 2 ( 2 ) 


A2(3) 

A2(4) 

A2(5) 

A2(6) 

A2(7) 

A2(8) 

A2(9) 

A2(l 0) 


the heat pipe characteristics array are: 

= Total heat pipe length (inches) 

= Number of heat pipe wall nodes (Floating 
Point Number) 

O 

= Heat pipe free volume to length ratio (in /in) 

= Reservoir volume (in ) 

= Gas NR value (Ft-lb/°R) 

= Reservior wick flag: 1.0 wicked reservoir 

-1.0 un-wicked reservoir 

= 1 .0 

= Reservoir heater power (Btu/Hr), -1.0 if unheated 

= Upper temperature setting (°F), (Heated 
reservoir system) 

= Lower temperature setting (°F), (Heated 
reservoir system) 
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A2(ll ) 
A2(12) 

A2(13) 

A2(14) 

A2(15) 

A2(l 6) 

A2(17) 

A2(18) 

A2(l 9) 


0.0 

Heat pipe identification number (Floating 
Point Number) 

A flag, 0.0 for steady state, 1.0 for 
transient solutions 

Working fluid freezing temperature 

A flag, 1.0 to activate ice plug option, 

0.0 otherwise 

Length of adiabatic section between evaporator 
and condenser 

Specified open artery (deprimed) power-length 
capac ( ty 

A flag, 1.0 to consider deprimed capacity, 

0.0 for unrestricted capacity 

Length of evaporator section. 


Sample Problem 

To illustrate the usage of subroutine VCHPDA, the steps 
involved in solving a simple problem are described in what 
follows. Consider the heat pipe radiator system sketched in 
Figure 3. The heat source is a 0.32 Kg aluminum block in which 
power is uniformly generated. The block is mounted on a 2,5 cm 
wide saddle that is attached to the heat pipe over a 30 cm-long 
section. The VCHP is a 1.25 cm O.D., 1,0 meter long aluminum/ 
ammonia heat pipe with a 0,1 litter gas reservoir. Short 
Al/SS/Al and Al/SS transition sections minimize the coupling 
by conduction between evaporator and condenser and between 
reservoir and condenser respectively. Power generated in the 
heater block is transported by the heat pipe to a radiator panel 
from which is radiated to a sink. 

The system is nodalized as shown in Figure 4. A listing 
of this model is shown in pages 45 through 47. The FORTRAN 
version of subroutine VCHPDA is listed in pages 48 through 67. 

Steady state, followed by transient solutions were obtained 
using SINDA. Selected solution outputs are shown in pages 68 
through 79. The results are summarized in Figure 5. 



Figure 3. Sample Problem Configuration 
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For the first hour after steady state conditions were 
reached, the VCHP system is shown to transfer 40 watts from 
the heat source to a sinK at 190°K. The vapor temperature 
remains constant at 296°K and the gas blocks about 3/5 of the 
condenser section. When the power is turned off, the vapor 
temperature can be seen to drop resulting in increase gas 
blockage of the pipe. At 1.5 hours the entire length of the 
pipe is gas blocked. Although the rate of heat leak from the 
evaporator is of only one percent of peak power (i.e., 0.40 
watt), the vapor, and consequently, the source temperatures 
undergo a substantial drop during the off period. This is 
the result of the rather small mass (0.32 Kg) of the source 
modelled in this problem. 

Continuous operation at low sink conditions, causes the 
inactive condenser to reach subfreezing temperatures and, as 
shown on the bottom of Figure 5, an ice plug forms at 2 hours, 
0.10 meters from the reservoir end of the pipe. The center 
sketch on the figure shows that a pressure differential de- 
velops across the ice plug. This barrier effectively decouples 
the evaporator from the gas reservoir and causes a large tem- 
perature overshoot when the power is turned on at 4.0 hours. 
Between this time and 5.0 hours, a large pressure differential 
is sustained across the ice plug. Following a step increase 
in sink temperature (5.0 hours) simulating a change in solar 
view factor, the ice plug thaws giving rise to the freezing 
blow by mechanism which is postulated to cause wicking failure. 
Continuous operation of the VCHP at 40 watts and deprimed ca- 
pacity (7.5 watt-meter) is seen to result in partial evaporator 
dry-out as the heatpipe adjusts its effective length to satisfy 
the imposed capacity constraint. 

References : 

Wanous, D. J., "Variable Conductance Heat Pipe Analysis Sub- 
routine, 1973 IR&D Project; Advanced Thermal Control", 

TRW IOC 8263.135-73-04. 
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00480 C /i?(15)=l. SHOUIL USED F'?P TRAMSTFNT RUNS ONLY, HOWEVPP 

00490 C 42(18)=!, COULL 3E USED FDR STEADY AS WELL AS FOR TRANSIENT 

00500 C PUNS. IN ADDITION A2(15)=l. AND A2(18)«l, COULD BE USED 

00510 C simultaneously FOR TRANSIENT PUNS. 


00520 C 
00530 C 

005A0 C I': THF PPEVni.'S VCH°2 IT WAS PEOUIRFO THAT A2(13) = .1 PF 

00550 C liSPL AS A OAHf-lNG IACTCr FOR STFAOY STATF RUNS AND THAT 
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00940 IF( A 2 ( 6 ).GT. 0 .) GO TO 102 

00950 CALL DIDFGKTWQ), Al.PVR) 

00960 C-n To 101 

00970 102 CALL OIPECK T?> Al, FVF ) 
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05120 GO TO 3000 
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APPENDIX A.b.l 

POTENTIAL FOR BUBBLE FORMATION IN CTS HEAT PIPES: 
FUNDAMENTALS 
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INTRODUCTION 

In the Anomalies Review and Program Plan of 19 October 1978, there 
is promised the following analysis: "Calculate number of critical 

size bubbles which can be generated due to supersaturation in rapid 
chill down. Consider both temperature and pressure effects". In 
what follows the basis for making such a calculation is reviewed. 

GAS SATURATION 

For dilute solutions of gas in liquid, Henry's Law can be invoked: 
The partial pressure gas i in the vapor phase P^. and the mole 
fraction X. or concentration c. of gas i in the liquid phase are 
linearly related 

^ " ^i h = [‘^i j c., i > 2 (1) 

In the vapor phase the total pressure P is given by Dalton's Law 
for an ideal gas mixture 

P - P, (2) 

t ■ 1 

where the vapor pressure of the solute species is given by 
Raoul t’s Law 

’’l ' ''„t <T) *1 ' C, (T) X, (3) 
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Often Henry's content (T) is so Urge that all X^'s 
other than X-j are very small, and X.j ■ 1 . The mole 
fractions in the liquid sum to unity, of course. 


i = 1 


i 1 


If the liquid contacts a gas mixture of specified total pressure P 
and specified ratios of noncondensible gas mole fraction, 


V. = P. 


T, ’’i - 


1 = 2 

First, one finds the mole fraction of noncnndensibl es in the 
liquid 


£ (Y,/C,) (P-C,) 

1 - 2 

n 

1 - X) 

i = 2 


The partial pressures in the vapor phase are 


= V, [ P - C, (1 . Xj) ■ 


The mole fractions in the liquid phase are 

*1 * ['■-Cl 0 -Xg)' 


CRITICAL BUBBLE SIZE 


The preceecHng relations are assumed to apply to gas at 
equilibrium within a critically-sized bubble when the total 
pressure in the bubble exceeds that in the surrounding 
liquid by an amount related to the surface tension a and bubble 
size r 

P - Pj, = L^JI) (9) 

Given a value of o (T), and r one can find P from 
Eq. (9) and then proceed as before to establish the set of 
mole fractions in the liquid immediately surrounding the 
bubble (the effect of X. on o is neglected). Conversely 
given a set of X^. , one can use Eqs. (1), (2), and (3) to 
find P and Eq. (9) to find r. 

PREVIOUS AND PRESENT STATE VARIABLES 

For convenience we choose to set the composition of a (bubbly) 
liquid by its '‘Previous state variables" T' and T' plus the Y. 
ratios. The quantity T'^ is the "Previous evaporator temperature". 
It sets the previous total pressure P' 

P' » p5,t‘T'e> 

The "Previous (gas-blocked) condenser temperature is T- ". It sets 
the previous vapor pressure according to Eq. (3) and, as explained 
below Eq. (4) in Eqs. (5) - (8), with P' and P'.j one can find 
the previous liquid composition X'^. 
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Present state variables are and T^. The quantity 
Tg is understood to fix pressure P^ by the relation 




The quantity T^ fixes a and P^ . 


( 11 ) 


CRITICAL BUBBLE SIZE IN AN INFINITE LIQUID 

Specification of the previous state variables, present state 
variables, and nonisothermal gas composition leads immediately 
to the critical size of a single bubble in contact with an 
infinite liquid. Variables T ', T ', and Y' . leads to a set of 

w W I 

X'^. which, in an infinite liquid, are identical to X^. . Then, 
as explained below Eq. (9), one can find P and r using the set 
of X^ and P^ (from T^) and a and P^ (from T^). 


cr 


" 2 o (T^) 


02 ) 


Where P^^ comes from Eq. (11) and P from Eq. (2) with P^ from 
Eqs. (1) and (3). 

NUMBER OF BUBBLES OF A SPECIFIED SIZE IN A FINITE LIQUID 

Let the mole fractions in one mole of liquid mixture be 
specified by the previous state variables. Imagine that 
the liquid is supersaturated, that is, there exists a finite 
r . Now choose a value r greater than r . Then Eq. (9), 

Cl C 1^ 

with P fixed by T and a by T , fixes total pressure P. Un- 
fortunately the gas composition in the bubble is unknown so 
that Eq. (7) cannot be directly applied. However, one can 
write that the moles of species i remaining in the liquid plus 
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those in the gas sum to the original amount of species i. 

(1 - N) + N «= X^.‘ (13) 


where N is the number of moles in the vapor phase. 
Introducing Eqs. (1), (2), and (3) gives 


X. + N P./P = X^.' 

X,. + (N/P) X. = X^ 

where C.j is understood to be Pgat^^c^' ^o^^ing for 
X^. gives 

X . = ^i 
^ 1 + (N/P) C. 


(14) 


(15) 


Multiplying both sides by C.(T )/P and summing over i 

I w 

gives a single equation which fixes the unknown N. 

s ’■r ( 16 ) 

P + N C, 

For a single-species noncondensible gas Eq. (16) may be 
solved explicity via the quadratic equiation. For a 
binary noncondensible gas mixture, one can employ a 
computer-aided binary search for the normalizing value 
of N. Once N is determined the number of bubbles follows 
immediately 

^ ^ u^c Bubbles (17) 

MfTTT 

where is the universal gas constant. 
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appendix a. 5. 2 

POTENTIAL FOR BUBBLE FORMATION IN CTS HEAT PIPES: 
SAMPLE CALCULATIONS 
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INTRODUCTION 


It is desired to calculate the number and size of critical size bubbles which 
can be generated due to supersaturation. In Part I the fundamentals were developed. 
Here sample calculations are made for two test cases (1) a condenser depressurization 
case where evaporator tenperature is dropped at constant gas -blocked condenser 
temperature and (2) a condenser chilldown where the gas -blocked condenser is cooled 
at constant total pressure. 

HENRY *S CONSTANTS 

Henr>''s constants are calculated from extrapolations of a solubility curve 
from Saaski. Saaski plots mole fraction of the gas in the liquid versus temperature 
on a log-log scale. Presumably the gas pressure is one atmosphere. Thus one has 
Xi for ?£ = 1» and since P^ = Cj^X^, C^ = P^/X^ = l/X^^. 

As a check, an Ostwald coefficient reported by Saaski for helium in methanol 
at 25°C is conpared to the graphed value. Ostwald coefficient o at 25°C is 
reported to be 0.036. The graph shows X^ = 6 x lO'^. Hence we expect from the 
graph that C^ = 1/6 x lO"^ = 16670 atm. Ostwald coefficient and Henr>^ number 
are related by 

_ 

^i a Ma 

c. . (0-785) (62.05) ( i29 8 ..1Jl . 

^ (32) (0.036) 

The check is very satisfactory. 

Table 1 gives Henry constant for helium and nitrogen at -100°C and -40°C. 

Hie values are hypothetical and extrapolated, hypothetical in that liquid is 
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freezes at a temperature of -98°C, and extrapoluted in that Saiaski's cxirves are 
extrapolated . The values are chosen merely to exenplify trends. 

TEST CASES < 

Table 2 shows the test-case conditions selected. Table 3 shows the property 
values needed for the test cases. 

PREV^IOUS STATE CO^IPOSITIONS 

As detailed in Part I the calculation commences with Eq.(I-6) and proceeds 
to Eq. (1-8). 


g 


[(72/^2^ ^ (y3/C3][P-q] 

i-[(y2./C2) + (y3/C3)]C;^ 


Xi=Cyi/Ci)[P-Ci(l-Xg) 

Table 4 gives the values. 
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TABLE 1 


Extraction of Henry's Constants From 
A Solubility Curve 

C = 1/x^ atm 


Gas 

Temperature 

Solubility 
in Methanol 

Henry’s Constant 
atm 

Helium 

-100C 

(Hypothetical) 

1 X 10'^ 

100000 


-40C 

(Extrapolated) 

2.6 X 10“^ 

38500 

Nitrogen 

-lOOC 

(Hypothetical) 

2,82 X 10"^ 

3550 


-40C 

(Extrapolated) 

2.69 X 10"^ 

3720 
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TABLE 2 


Ffypothetical Test Cases 


Test Case 1 Depressurization from Drop in Evaporator Temperature 
Previous State Variables 

Evaporator Temperature, T ' 

w 

Gas Blocked Condenser Tenq^eratiore, T^’ 


Noncondensible Gas Composition, y 


n,i 


49C(120F) 

-40C(-40F) 

90% N 2 - 10% He 


Present State Variables 

Evaporator Temperature, T 

c 

Gas Blocked Condenser Temperature, T 


21C(70F) 

-40C(-40F) 


Test Case 2 ChiHdown in Gas -Blocked Condenser 
Previous State Variables 

Evaporator Temperature, T ' 21C(70F) 

Gas Blocked Condenser Temperature, T^' -40C(-40F) 

Noncondensible Gas Composition, ^ 90% N 2 - 10% He 

Present State Variables 

Ev'aporator Temperature, T^ 

Gas Blocked Condenser Tenperature, T^ 


21C(70F) 

-100C(-148F) 


TABLE 3 


Vapor Pressure and Surface Tension 


T 

C 

^sat“^l 

psia 

T 

C 

lbf/£t 

49 

7.62 

-40 

2.1 X 10"^ 

21 

1.96 

-100 

2.4 X 10'^ 

-40 

0.02870 



-100 

0.00012 
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Previous State Conditions 
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Critical Bubble Size in an Infinite Liquid 


As explained in Part I, to coupute the ciritical bubble size, one conputes 


P 





and then 


cr 


2a (Tg) 


Table 5 shown the results. For depressurization (case 1) the critical bubble 

is small, 1.6 ym in radius, the size of nuclei expected to be active in nucleate 
surface boiling. However, the gas conposition near the surface exposed to the 
gas wrould be more nearly in equilibrium. That is contact vdth the wick might 
be close to the ccmposition assumed. For chilldown (case 2) the increase in 
Henry’s constant for helium results in supersaturation, but the potential for 
nucleation is much less. A fifty -micron-radius nucleus would be required for 
bubble formation. 


Number of Bubbles of a Specified Size 

From Part I the equation which governs the amount of gas N within bubbles 
is 


1 


« 


Cl Xj^ 

P+FiC^ 


f 


2 ^2 
P+5^ 


^3 ’‘S 


vdiere the total pressure within the bubble is 


P 


' '’sattV- * 


2a 

r 
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For case 1 and r ■ 50.8 ym the following values pertain: 


P » 1.96 + = 2.135 psia « 0.1453 atm 

50.8/25400 

CjX^’ «= (0.001963) Cl. 000) «= 0.001953 atm 
^2^2* “ (3720) (1.251 x lO"^) = 0.465 atm 
C 3 X 3 ' = (38500) (1.34 X 10"^) * 0.052 atm 
N = 8.96 X 10 ’^ moles of bubbles/mole of original liquid 
The number of bubbles per mole is then 

w = N RiiTr ^ (8.96 X 10 '^) (82.05) (233.15) 

^ P(| (0.1453) (|^ tt) (50.8 X lO'"^)^ 

= 2.15 X lo”^ Bubbles/g-mole of liquid 

Taking to be approximately 0.8 [g/cm^3/32[g/g-mole] gives 
= 537000 Bubbles/cm^ of liquid 

Conclusions 

The sample calculations show the possibility for bubble formation upon 
( 1 ) depressurization by reduction in heat pipe loading and reservoir chilling 
and (2) condenser chilldown. Depressurization shows greater potential, for 
537000 bubbles/cm3 are found possible, compared to only one of equal size 
from condenser chilldown. 
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APPENDIX A. 5. 3 

BUBBLE NUCLEATION EXPERIMENTS 
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' As part of the experimental program to investigate the potential of 

I various medianisms to cause artery depriming, a series of experiments have * 

been conducted to examine bubble nucleation in methanol. Tlie experiments | 

were to deteimiine v\rhether gas bubbles could be generated in the bulk of ) 

liquid methanol saturated with either helium or nitrogen gas as it undergoes I 

tojnperature and/or pressure reduction. 

Prior to the experiments, a theoretical analysis of the potential for 
bubble formation in Cl’S Jicat pipes had sliown that large numbers of bubbles 
could be generated in methanol due to supersaturation resultirig from temper- 

• ature and pressure reduction under conditions similar to those prevailing 

* prior to the anomalies. An objective of these experiments was to verify, ■* 

' at least qualitatively, tlie theoretical results. 

i In addition the ejqseriments considered the potential of a mesh screen 

1 artery to provide bubble nucleating sites. , 

i The experimental set-up consisted of a glass flask half filled with 

50cc of spectral grade metlianol and instrumented to allow continuous 
monitoring of temperature and pressure. A sketdi of the apparatus is 
attached. A needle valve located between the flask and a vacuum pump 
' was used to control the pressure level or tlie rate of pressure reduction 

j of the liquid. Cooling of the liquid was attained by immersion of the 

test flask in liquid nitrogen. This teclmique, however, did not allow ^ 

for arbitrary control of the cooling rate. The liquid was saturated by 

bul)blirig either nitrogen or helium gas through the liquid using a frit 

glass tube. Tliis process was allowed to be continued for at least two 

hours . . 

A t>T>ical pressure reduction experimental sequence was as follows: 

1. Saturate mctJianol with either niti'ogcn or lielium at ambient 
conditions. 

2. Reduce temperature at atmospheric pressure. Two temperature 
levels, 21°C and -40‘^C, were used. 

3. Drop into the liquid a dry section of mes]\ screen artery. 
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4. Reduce pressure from 14.7 psia to 4 psia. In some cases this 
pressure drop was accomplished in 10 seconds and in others in 
5 minutes. 


A typical temperature reduction experimental sequence was as follows: 

1. Saturate methanol at ambient conditions. 

2. Drop into the liquid a dry section of mesh screen artery. 

3. Reduce pressure at ambient temperature. Pressure levels of 
14.7 psia, 10 psia, and 4 psia were used. 

4. Reduce temperature. In some cases the temperature was reduced 
from 21°C to -40°C in approximately 20 minutes. This was accom- 
plished by placing the flask near the liquid nitrogen surface 
contained in a Dewar flask where cooling of the methanol occurred 
by natural convection. In other cases, the liquid was rapidly 
chilled by immersing the flask in liquid nitrr'gen for approximately 
10 seconds at which point the liquid on the bottom of the test flask 
and inside the artery froze. 


The experiments were repeated several times and the results were found 
reproducible. No bubbles were observed as the liquid, originally at a set 
temperature level, underwent pressure reduction. Sunilar results were obtained 
from temperature reduction experiments provided the liquid te-mperature did not 
drop below the freezing point (-9S°C) . 

Experiments in which the liquid partially froze however, yielded significantly 
different results. A large niomber of small bubbles were observed streaming from the 
surface of the tlmwing ice. As the melting process went to competion tlic bubbles 
were reduced to a few originating from the bottom surface of the flask and from the 
outter surface of the mesh screen artery. About a dozen small bubbles were observed 
trapped inside the arteries after the ice melted. These bubbles were observed later 
to coalesce forming fewer but larger bubbles which continue to grow. The ultimate 
size of these bubbles depended on the pressure of the system. For example, at 4 psia 
the remaining bubbles inside the arter)' continue to grow for approxinatcly 10 minutes 
as tlie liquid tempei'ature rose from -40°C to 0°C, ac which point the size of the 
bubbles was such that essentially all the liquid in the artery was displaced. 

As the re.sult of those experiments, freezing of supersaturated methanol in the 
arteries has been identified as a potential mechanism for bubbl.e formation in she 
arteries owing to tlie fact that the ice surface is an excellent source of nucleating 
sites. 
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APPENDIX A.5.4 

GLASS HEAT PIPE BUBBLE NUCLEATION/MIGRATION EXPERIMENTS 


A series of experiments were performed with an existing 1.07 meter 
long glass heat pipe. A cross section of this heat pipe is shown in Fig- 
ure 1. The pipe contains a slab wick with a CTS-type artery attached to 
one side. A heater and cooling loop are attached to the other side of the 
wick in the evaporator and condenser sections, respectively. This arrange- 
ment permits observations on the behavior of the artery in an operational 
heat pipe. 

The heat pipe was gas loaded with a 90% nitrogen — 10% helium mixture 
at a pressure equivalent to the conditions in the CTS pipes. The experi- 
ments simply involved visually observing the artery behavior as a result 
of freezing the methanol within it by passing liquid nitrogen through the 
cooling loop, and subsequently thawing the methanol by terminating the LN 2 
flow. 


The experiment was repeated several times with the following results: 

• The methanol froze to an opaque white solid (frost), indi- 
cating it was full of gas bubbles. 

• As the methanol thawed and warmed, numerous (up to 33) gas 
bubbles of varying size were observed along a IZ^inch length 
of the artery. 

t The bubbles slowly shrank and ultimately disappeared as the 
gas within them redirsolved into the liquid and diffused 
through the artery wall into the surroundir-q vapor core 
under the pressure gradient established by their curvature 
and surface tension. 

• The time required for the bubbles to disappear varied enor- 
mously, depending on their size. Very small bubbles 

(d '^'l/4 artery dia) dissolved within a few hours, while a 
large sausage-shaped bubble which filled the artery (2 
artery diameters long) required several weeks. 

It was believed that incorporation of these ice-generated bubbles into 
the active condenser could cause artery derriming. Accordingly, additional 
tests' were conducted and a number of deprimings were observed, but they 
were sporadic (probabilistic in nature). On some occasions, as the active 
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condensation front was brought into the previously gas-blocked condenser, 
no depriming occurred. On two occasions, freezing generated a bubble in 
the adiabatic section; on both these occasions, advancement of the front 
to the bubble caused depriming. On another occasion, a rather large bubble 
formed in the center of the condenser, and advancement of the front to this 
bubble caused depriming. 

These results conclusively demonstrated that: 1) control gas is 

liberated from saturated methanol every time the heat pipe goes through 
a freeze-thaw cycle, and 2) the number, size, and durability of gas bubbles 
generated within the arteries is a statistical phenomena influenced by bub- 
ble coalescence. Arterial bubbles are known to lead to depriming if they 
are convected into the active region of the pipe under high load conditions. 
Thus, it if. clear that freezing and thawing of the condenser can, but does 
not necessarily, lead to artery depriming, depending on subsequent history. 
That is, do the bubbles dissolve and diffuse away before they are con- 
vected into a high stress region where they would deprime the artery? 

This mechanism appears to be a prime candidate for explaining the CTS 
anomalies. It is consistent with their seasonal occurrence (eclipse sea- 
sons may be necessary to experience condenser freezing), sporadic occur- 
rence (statistical nature of bubble population) and the lack of freezing 
prior to the anomaly on day 253 (bubbles were generated during day 252). 
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APPENDIX B.1 


SPHERICAL BUBBLE MODEL COMPUTER PROGRAM LISTING 


This appendix contains the main program and sample input 
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05560 FMAX » -BIG 

05570 PMAX » -BIG 

05560 NP = 0 

05590 return 

05600 END 
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Sample Input 


■4.BUBDAT 

C0S11»2.5091E-12 
C0S21=2«9665 
C0S12=2,6367E-^ 
COS22»0.0 
Df'U»5.92E-3 
=<«1,03A ,2 
DF12‘^?-.09E-3 
DF22--1070. 5 
NG1 = 1 
NG2«2 

NSTPMX=8000 

IDPLOT-30 

R10»0.02 

R20=0.08 

NDIV=90 

GAMA--2 .0 

DT»5 . 

IDPRNT-200 

RFILLET=10. 

FBYl-1.0 

FBY2»0. 

FEY1=0.1 

FEY2*0.9 

EMW«32.0 

TCN0T=2&0. 

TEN0T=300. 

RATE=0.0 

$6ND 
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APPENDIX B.2 

SN009 HEAT PIPE CYCLIC TEST DATA 

This appendix presents plots of temperatures along the SN009 heat pipe 
versus time for selected freeze/thaw test cycles not shown in Section 3.3. 
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TEMPERATURE SENSOR LAYOUT 



TIME (MINUTES) 

Figure B.2-1 Temperature Histories During SN009 Heat Pipe Cyclic Test No 




EMPERATURE SENSOR LAYOUT 



gure B.2-2 Temperature Histories During SN009 Heat Pipe Cyclic Test No. 







TIME (MINUTES 


TEMPERATURE SENSOR LAYOUT 



TIME (MINUTES) 

Figure B.2-5 Temperature Histories During SN009 Heat Pipe Cyclic Test No. 10 





TEMPERATURE SENSOR LAYOUT 



TIME (MINUTES) 

tojiCS During SN009 Heat Pipe Cyclic Test Ho. 14 



TEMPERATURE SENSOR LAYOUT 



TIME (MINUTES) 

Figure B.2-8 Temperature Histories During SN009 Heat Pipe Cyclic Test No. IfJ 






TEMPERATURE SENSOR LAYOUT 



Figure B.2-9 Temperature Histories During SN009 Heat Pipe Cyclic Test No. 17 
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Figure B.2-11 Temperature Histories During SN0Q9 Heat Pipe Cyclic Test No. 25 



pe Cyclic TPSt Nn. 



